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ABSTRACT 


The  Ballistic  parameters  were  determined  analytically  and  experimentally 
for  the  perforation  of  aluminum  alloy  and  titanium  alloy  plates  by  normal  im¬ 
pact  of  steel  cylinders.  Analytical  predictions  of  the  perforation  parameters 
were  prepared  on  the  basis  of  theoretical  models  of  the  perforation  process. 
Experimental  determination  of  the  critical  perforation  and  residual  velocities 
discriminated  among  the  various  models  and  further  experimentation  verified  the 
analytical  description  of  the  perforation  process. 

One  model  of  perforation  was  found  which  adequately  describes  the  ballistic 
parameters  for  the  normal  impact  of  steel  cylinders  into  aluminum  alloy  and 
titanium  alloy  plates  for  impact  velocities  up  to  5000  ft/sec  and  plate  thick¬ 
nesses  up  to  0.25  inches. 
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SECTION  1 


INTRODUCTION 

1.1  Objectives 

The  results  of  an  analytical  and  experimental  investigation  of  the  per¬ 
foration  of  thin  metallic  plates  by  the  normal  impact  of  steel  cylinders  are 
presented  in  this  report.  This  study  follows  a  previous  research  program  which 
investigated  the  basic  physical  mechanisms  of  deformation  and  failure  of  metal¬ 
lic  plates  under  ballistic  impact. 

Two  previous  reports  described  the  results  of  analytical  and  theoretical 
investigations  into  the  nature  of  the  perforation  process.  The  subject  matter 
was  highly  mathematical  and  concerned  itself  mostly  with  descriptions  or 
attempts  at  descriptions  of  the  Physical  phenomena  that  occur  during  impact. 
Analytical  and  numerical  solutions  for  the  propagation  of  stresses  into  single 
and  composite  plate  systems  were  given  by  the  linear  theory  of  elasticity; 
for  a  simplified  geometry  of  the  single  plate,  a  solution  to  the  propagation 
of  stresses  in  an  elastic-plastic  plate  was  obtained.  Supplementary  investi¬ 
gations  of  the  properties  of  materials  needed  for  evaluation  of  fracture  and 
plastic  flow  were  also  performed. 

The  previous  analytical  work  delineated  several  mechanisms  through  which 
the  plate  may  fail  under  ballistic  impact.  Each  mechanism  of  plate  failure 
defines  a  model  of  perforation.  Each  model  of  perforation  defines  a  character¬ 
istic  critical  perforation  velocity  and  characteristic  residual  velocities. 

These  ballistic  parameters  are  functions  of  plate  geometry,  projectile  geometry, 
projectile  material  properties  and  plate  material  properties. 


The  objective  of  this  studywas  to  determine  which  of  the  models  of  per¬ 


foration  is  descriptive  of  the  perforation  of  thin  metallic  plates. 

1.2  Method  of  Approach 

To  accomplish  the  objective  of  this  study,  the  various  modes  of  failure 
were  delineated.  The  critical  perforation  velocity  and  residual  velocities 
were  predicted  for  the  perforation  of  single  metal  plates  under  normal  impact 
by  right  circular  cylinders.  In  addition,  an  experimental  determination  of 
these  quantities  vf  as  made. 

The  investigation  of  the  perforation  parameters  was  composed  of  two  phases 

Phase  I:  Critical  velocity  predictions  were  prepared  for  all  the 

models  of  perforation  for  a  selected  cylindrical  pro¬ 
jectile  and  two  different  plate  materials.  Predictions 
were  made  for  several  possible  models  of  perforation. 

The  only  parameter  allowed  to  vary  was  the  plate  thick¬ 
ness.  Tests  were  then  conducted  to  experimentally 
determine  the  critical  perforation  velocities  under 
these  conditions.  Comparison  of  the  analytic  predictions 
and  the  experimental  results  showed  that  two  models  of 
perforation  used  in  comhinatLon  fitted  the  experimental 
data.  Another  model  of  perforation  also  fitted  the 
experimental  data.  Each  of  the  models  has  distinguishing 
characteristics  which  suggested  the  subsequent  choice  of 
the  projectile  and  plate  parameters  for  the  second  phase 
of  the  investigation. 
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Phase  II:  Further  experimental  determinations  of  the  critical 

perforation  velocities  for  one  of  the  plate  materials 
were  made,  using  cylindrical  projectiles  of  different 
sizes.  The  projectiles  used  in  these  tests  were  chosen 
using  the  experimental  data  and  the  theoretical  models  of 
perforation.  The  projectiles  were  selected  to  provide  a 
distinguishing  test  among  the  perforation  models  that 
might  adequately  explain  and  describe  the  process. 

1.3  Summary  and  Conclusions 

The  analytical  and  experimental  portions  of  the  research  program  were 
carried  out  as  indicated  in  Section  1.2.  The  previous  analytical  work  was 
reviewed  and  several  possible  modes  of  failure  of  the  plate  under  ballistic 
attack  were  delineated.  Analytical  expressions  were  derived  for  the  critical 
perforation  velocity  and  the  residual  velocities  for  each  of  the  several  modes 
of  failure  of  the  plate  system.  The  ballistic  parameters,  in  particular  the 
critical  perforation  velocity,  were  computed  for  the  normal  impact  of  right 
circular  cylinders  of  hard  steel  into  two  different  metallic  alloys.  The 
critical  perforation  velocity  and  the  residual  velocity  of  the  projectile 
as  a  function  of  the  striking  velocity  of  the  projectile  were  experimentally 
determined  for  the  impact  of  steel  cylinders  into  two  metallic  alloy  plates. 

Experimental  data  were  compared  with  theoretical  predictions  in  both 
phases  of  the  work.  The  ’’modified  shear  plug"  mode  of  failure  of  the  plate 
was  deemed  to  be  the  most  descriptive  for  prediction  of  the  critical  per¬ 
foration  velocity  for  the  range  of  impact  velocities,  plate  thicknesses,  plate 
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materials  and  projectiles  tested.  This  model  gave  accurate  predictions  for  the 
critical  perforation  velocities.  The  residual  velocity  versus  striking  velocity 
curves  as  predicted  by  the  model  are  in  good  agreement  with  the  experimental 
observations.  Further,  the  model  is  readily  adaptable  to  the  computation  of  the 
perforation  characteristics  of  projectiles  of  different  nose  shapes. 

Under  slightly  different  conditions  of  impact  of  the  projectiles,  and  at 
higher  impact  velocities,  the  mode  of  the  failure  of  the  plate  may  change. 

This  point  will  be  discussed  in  the  companion  volume  to  this  report  in  conjunction 
with  the  prediction  of  the  critical  perforation  velocities  of  military  projectiles. 

1.4  Organization  of  the  Remainder  of  the  Report 

In  Chapter  2, a  brief  non-mathematical  summary  of  the  physics  of  the  per¬ 
foration  and  deformation  of  a  thin  metallic  plate  under  ballistic  impact  is 
given,  together  with  a  list  of  the  critical  velocity  formulae  for  the  several 
modes  of  failure. 

Chapter  3  presents  the  comparison  between  theory  and  experiment  for  the  per¬ 
foration  of  the  plate  by  normal  impact  of  right  circular  cylinders. 

Two  appendices  are  attached  to  the  report.  The  first  appendix  describes 
the  experiments.  The  second  appendix  presents  a  detailed  description  of  the 
physics  of  impact  and  perforation,  the  mathematical  models  used  for  the 
description,  the  detailed  findings  of  the  analytical  research  conducted  on  this 
problem,  and  the  complete  derivation  of  the  model  for  perforation  which 
best  fits  the  experimental  data. 
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A  companion  volume  contains  predictions  of  the  critical  and  residual 
velocities  expected  for  normal  impact  of  various  steel  projectiles  into  several 
plate  materials.  The  projectiles  are  fragment  simulators  and  sharp  steel  bullets 
of  various  dimensions.  Steel,  titanium  alloy  and  aluminum  alloy  plates  are  the 
plates  considered.  Description  of  the  application  of  the  theory  to  the  determi¬ 
nation  of  an  optimum  armor  system  is  also  given.  Comparison  of  the  predictions 
with  experimental  data  on  the  perforation  of  the  metal  plates  by  the  projectiles 
is  given. 


t; 


SECTION  2 


THE  PHYSICAL  PHENOMENA  OF  THE  PERFORATION  PROCESS 

This  section  summarizes  the  basic  sequence  of  events  in  the  plate  during 
impact  and  perforation,  delineates  the  models  of  fracture  to  which  the  plate 
is  subjected,  and  enumerates  the  perforation  predictions  for  each  mode  of 
failure.  Appendix  B  contains  a  detailed  and,  for  the  most  part,  non- 
mathematical  description  of  the  physics  of  perforation. 

Consider  a  right  circular  cylindrical  projectile  of  radius  a,  and  length  L 

approaching  a  plate  of  uniform  thickness,  h.  The  projectile  is  oriented  so  that 

its  axis  is  perpendicular  to  the  plate.  The  projectile  approaches  the  plate 

with  a  velocity  V  ,  which  is  perpendicular  to  the  surface  of  the  plate.  There 

s 

is  no  lateral  velocity  component.  When  the  projectile  strikes  the  plate,  the 
entire  flat  end  of  the  cylinder  will  strike  the  flat  surface  of  the  plate. 

The  density  of  the  plate  is  and  the  density  of  the  projectile  is  PQ. 

(See  Figure  1.) 

The  physical  processes  that  occur  immediately  after  the  projectile  impacts 
the  plate  are: 

1.  A  large  stress  is  generated  by  the  impact  at  the  interface 
between  the  projectile  and  the  plate.  This  stress  is,  to 
first  approximation,  proportional  to  both  the  acoustic  impedance 
of  the  plate-projectile  combination  and  the  impact  velocity. 

The  impact  stress  is  a  monotonic  function  of  the  impact  velocity. 
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Figure  1  GEOMETRY  OF  IMPACT  AND  DEFINITION  OF  TERMS 


2.  Two  stress  waves  are  propagated  into  the  plate:  a  fast 
dilatational.  wave  and  a  slower  shear  wave.  The  stress  in 
the  first  wave  is  compressive.  The  stress  in  this  first  wave 
rises  very  rapidly  and  then  decays  very  rapidly.  Immediately 
behind  this  spike  a  large  stress  is  left,  but  not  as  large  as 
the  stress  in  the  spike  The  second  stress  wave  is  the  shear 
wave.  For  the  normal  impa-  of  a  cylinder  into  a  half- space, 
the  transient  stresses  in  this  wave  are  small.  After  the 
passage  of  the  shear  front  in  a  half-space,  the  quasi-static 
or  steady  state  solution  is  established  very  rapidly. 

3.  These  stress  waves  reflect  from  the  back  surface  of  the 
plate,  which  is  a  free  surface.  The  first  of  these 
reflected  stress  waves  is  tensile.  Meanwhile,  two  stress 
waves  are  propagating  into  the  projectile  and,  after  trav¬ 
ersing  the  projectile,  will  reflect  from  the  end  and  return 
to  the  impact  surface.  These  stress  waves  will  reflect  from 
the  original  interface,  the  reflected  stress  wave  again  trav¬ 
ersing  the  plate.  These  reflections  will  continue  for  a 
long  time. 

U.  After  the  stress  wave  has  completed  the  first  full  traverse 
of  the  plate,  the  normal  stress  becomes  small.  The  portion 
of  the  plate  under  the  projectile  is  in  motion  with  a  normal 
velocity  slightly  less  than  the  impact  velocity.  High  radial 
and  tangential  stresses  are  developed  in  the  plate.  High  shear 
stresses  are  generated  in  the  plate  near  the  lip  of  the  projectile. 


8 


In  response  to  the  stresses  that  are  generated  in  the  plate  system,  the 
plate  may  possibly  fail  by  one  cf  several  mechanisms.  They  are: 

1.  Fracture  in  the  Direct  Wave  -  The  plate  may  fracture  in  response 
to  the  high  compressive  stresses  that  are  generated  in  the  stress 
spike  near  the  first  stress  wave  front  that  propagates  into  the  . 
plate . 

2.  Radial  Fracture  in  the  First  Stress  Wave  -  After  the  stress  spike 
near  the  wave  front  of  the  longitudinal  wave  has  passed,  the  radial 
stress  in  the  plate  becomes  tensile  while  the  normal  stress  remains 
compressive. 

3.  Spall  -  The  plate  may  fracture  in  response  to  high  tensile 
vertical  stress  generated  by  the  reflection  of  the  stress  wave 
from  the  back  surface. 

4.  retailing  -  After  the  passage  of  the  first  stress  waves,  the 
plate  may  fracture  in  response  to  high  radial  and  circumferential 
tensile  stresses  that  are  developed  in  the  plate  near  the  lip  of 
the  projectile. 

5.  Plugging  -  After  the  passage  of  the  first,  stress  wave,  the  plate 
may  fail  in  response  to  the  large  shear  stress  that  is  developed 
near  the  edge  of  the  projectile. 

The  determination  of  the  ballistic  performance  cf  a  plate  or  .aminated 
plate  system  under  impact  by  a  projectile  has  two  basic  problems. 

1.  Does  the  plate  or  laminated  plate  system  fail  or  break? 

Does  the  projectile  stop  or,  at  least,  how  much  does  it  slow  down? 

o 


2. 


To  answer  the  first  question,  the  transient  stress  pattern  in  the  plate  that 
arises  from  the  sudden  generation  of  stress  at  the  contact  surface  must  be  known, 
together  with  the  criterion  for  failure  for  the  material  in  response  to  these 
time-dependent  stress  fields. 

To  answer  the  second  question,  we  need  to  know  the  magnitude  of  the  forces 
resisting  the  motion  of  the  projectile  and  how  long  they  are  maintained.  The 
magnitude  of  the  resisting  force  is  given  by  the  impact  stress  generated  over 
the  contact  area  of  the  plate  and  projectile.  This  force  will  be  maintained 
until  the  plate  breaks  or  until  modified  by  the  multiple  reflections  of  stress 
waves  that  occur  in  a  plate  system. 

The  critical  velocity  for  perforation  is  given  by  answering  the  second 
question,  while  the  critical  velocity  for  spalling  without  perforation  is 
given  by  the  answer  to  the  first  question.  In  the  following  paragraphs,  the 
perforation  limits  and  the  residual  velocity  of  the  projectile  will  be  given 
as  a  function  of  natural  parameters  for  the  various  modes  of  failure.  These 
critical  perforation  velocities  and  residual  velocities  were  calculated  on  the 
basis  of  stress  patterns  developed  in  the  plate.  For  the  first  three  modes  of 
failure  the  stress  wave  propagation  in  the  plate  was  assumed  linear  elastic. 
Further,  stresses  were  computed  on  the  assumption  that  the  projectile  did  not 
deform.  Shatter  of  the  projectile  is  ignored  here. 

The  material  parameters  and  their  effect  on  the  various  modes  of  failure 
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are  listed  below: 


1.  Fracture  in  the  Direct  Wave  -  The  material  parameters  involved 


in  the  direct  wave  fracture  at  the  back  side  of  the  plate  are  the  Brinell 
hardness  of  the  plate  (or^  equivalently,  the  yield  strength),  the  density  of 
the  plate,  and  the  compressive  ultimate  stress.  Increases  in  any  of  these 
reduce  the  penetration  due  to  fracture  in  the  direct  wave.  The  critical 
velocity  for  perforation  is 


where 


where 


V 

c 


1  -  2V 
1  -  V 


,1/2 

Z 


0,  is  the  tensile  stress 
ty 

Z  is  the  acoustic  impedance,  given  by 


Z 


P1C1P0C0 

pici+poco 


V 

s 


Pq,P^  are  the  densities  of  projectile  and  plate,  respectively. 

CQ,C  are  the  longitudinal  sonic  velocities  of  projectile  and 
plate,  respectively. 

2  2 

f  is  a  geometric  factor  given  by  i  +  h  /a 

h  is  the  plate  thickness 

a  is  the  projectile  radius 

V  is  the  striking  velocity 
s 

vis  Poisson’s  ratio  of  the  material 


2.  Radial  Fracture  -  The  critical  velocity  for  failure  in  •;he  reflected 
stress  wave  may  occur  when  the  radial  stress  exceeds  the  yield  stress.  For  a 
metal  plate  under  normal  impact  from  a  rigid  circular  cylinder,  this  critical 
velocity  is  given  approximately  by 


V 

c 


Z 


g(l  -  v  )f 

(i  _  v )  +  tv:1/2 


3.  Spall  -  The  material  and  geometrical  parameters  that  govern  the  spall 
are  the  thickness  to  bullet  radius  ratio  and  the  tensile  strength  of  the  plate 
in  the  direction  of  the  projectile  approach.  Tne  critical  perforation  velocity 
is  proportional  to  the  tensile  yield  strength  and  a  geometric  function.  The 
critical  velocity  for  failure  by  spall  is  determined  by  the  vertical  or  normal 
stress  and  is  bounded  and  closely  approximated  by  the  expression 


V 

c 


If  the  time  to  fracture,  t.^,  at  a  given  stress  level  is  taken  into  account 
in  this  mode  of  failure,  the  impact  stress  required  to  cause  failure  is  raised. 
The  critical  perforation  velocity  is 


with 


V 


c 


1  -  2y 
1  -  v 


10Vi 
hz  fi/2  e  ~ 

Z 


=  time  to  fracture 
=  longitudinal  sonic  velocity 


Computations  using  t^,  =  10  sec  were  used  at  first,  which  agrees  with  the 
experimental  data  for  time  to  fracture  in  terms  of  pure  aluminum.  The  effect  of 
time  delay  to  fracture  does  not  significantly  alter  the  predictions  of  the 
critical  perforation  velocity. 


The  stresses  in  the  spike  decay  extremely  rapidly  by  dissipation  in  plastic 
flow.  Explicit  expressions  for  the  decay  of  a  stress  wave  by  plastic  flow  have 
been  computed  for  the  propagation  of  a  spherical  wave  (Fugelso  1962)  and  a 
cylindrically  expanding  wave  (Fugelso,  Arentz,  and  Davidson  1962).  For  impact 


stresses  far  above  the  yield  stress  for  the  distance  the  stress  wave  must  propa¬ 
gate  in  the  present  problem,  the  stresses  should  decay  by  a  factor  of  1.5  to  2 
by  the  time  the  stress  wave  reaches  the  back  of  the  plate .  Thus  the  critical 
perforation  velocity  estimate  should  be  raised  by  a  factor  of  1.5  to  2. 

4.  &  5.  Plugging  and  Petalling  -  If  the  plate  does  not  fail  by  fracture 
as  the  first-stress  wave  passes,  then  the  failure  of  the  plate  by  plugging  and 
petalling  must  be  considered.  Plugging  is  failure  of  the  plate  by  shear 
under  the  rim  of  the  bullet,  whereby  a  plug  of  plate  material  is  ejected. 
Petalling  occurs  when  the  plate  fails  in  response  to  the  tensile  membrane 
stresses  that  are  generated.  The  critical  velocities  for  perforation  are 
similar  for  the  petalling  and  plugging  failure  modes.  The  static  components 
of  the  resisting  force  to  the  projectile  motion  axt  given  by  elementary  plate 
theory.  In  the  case  of  the  shear  fracture  (plugging  model)  the  resisting 
force  of  the  plate  was  taken  as  the  maximum  shear  stress  developed  in  plastic 
flow.  This  is  integrated  across  the  plate. 


a 


h 

°ty 


The  resisting  force  to  the  projectile  motion  is  given  by  the  maximum  stress 
resultant  in  the  plate. 


o  = 


with 


tensile  yield  stress  of  the  plate  material 


The  critical  perforation  velocity  here  is  given: 


V  =  Xl/2(eY  -  l)1/2 


for  plugging 

for  petalling 

For  both  mechanisms, the  critical  perforation  velocity  is  a  monotonic 
function  of  the  tensile  yield  stress  of  the  plate  and  the  plate  thickness  to 
projectile  length  ratio. 

The  mechanisms  are  normally  restricted  from  thin  to  moderately  thick 
elastic  plates  (h/a  <  l)  and  become  the  dominant  mechanisms  of  failure  for 
thin  plate  and  low  impact  velocities. 

One  further  revision  of  the  shear  model  was  made  and  calculations  per¬ 
formed.  The  static  component  of  the  resisting  force  of  the  plate  was  taken 
as  O  (the  tensile  ultimate  strength)  times  the  impacted  area.  The  resisting 
force  was  chosen  to  be  independent  of  plate  thickness  (unlike  the  previous 
model,  which  had  the  static  resisting  forces  proportional  to  the  plate  thick¬ 
ness).  The  resulting  formula  for  the  critical  velocity  for  the  perforation 
of  a  thin  plate  under  normal  impact  of  a  right  circular  cylinder  is 


2  P±h 

a,  , 

JSL 

-I 

p  l  1 

1  a  1 

Urr 

h 

Pi 

a 

V 

c 

X 

Y 


2Plh 


PfL 


'4 


where  0  is  the  tensile  ultimate  stress 
tu 

P(yP1  are  the  densities  of  the  projectile  and  plate 
h  is  the  thickness  of  the  plate 
L  is  the  length  of  the  projectile 

A  complete  derivation  of  the  "modified  shear  plug"  model  is  given  in  full 
detail  in  Appendix  B. 

Figures  2  and  3  show  the  functional  behavior  of  each  of  the  models  of 
perforation  as  a  function  of  increasing  (h/a)  ratios.  From  the  inspection 
of  Figure  2,  spall  can  immediately  be  eliminated  as  descriptive  of  the  critical 
perforation  velocity.  For  thin  plates  and  sufficiently  long  projectiles,  the 
c:  '.tical  velocity  for  spall  failure  has  a  minimum  near  —  =  1.5  and  increases 
without  bound  for  both  thicker  and  thinner  plates.  The  radial  fracture  and 
fracture  in  the  direct  wave  model  are  almost  identical  in  their  predictions 
of  the  critical  perforation  velocity.  Figure  3  shows  typical  behavior  of  the 
critical  perforation  velocity  for  the  petalling,  plugging  and  modified  shear 
model  as  functions  of  the  plate  thickness  to  projectile  radius  ratio.  The 
graph  shows  the  predicted  critical  perforation  velocity  for  7075- T6  aluminum 
plate  being  impacted  by  a  .22  caliber  4340  steel  cylinder. 
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(Plate  Thickness )/ (Projectile  Radius)  b/a 


Figure  2  DEPENDENCE  07  CRITICAL  PERFORATION  VELOCITY 
ON  PLATE  THICKNESS 


1.6 


Critical  Perforation  Velocity  (Ft/Sec) 


SECTION  3 


CALCULATIONS  AND  EXPERIMENTAL  RESULTS 


3.1  Experimental 

Predictions  of  the  critical  perforation  velocities  for  the  various  models 
of  perforation  were  made  for  two  plate  materials  under  normal  impact  by  right 
circular  cylinders.  The  plate  materials  were  the  titanium  alloy  Ti  5  A1  2.  5  Sn 
and  the  aluminum  alloy  7075-T6.  The  impacting  cylinders  were  4340  steel 
(normalized) . 

Experiments  were  performed  to  obtain  residual  velocity  vs.  striking  velocity 
curves  and  critical  perforation  velocities.  These  results  were  used  to  establish 
and  to  distinguish  the  various  perforation  models.  Three  types  of  projectiles 
were  shot  at  various  thicknesses  of  two  types  of  target  plates.  The  projectiles 
and  the  plates  used  in  the  experimental  phases  of  this  study  are  described  in 
Tables  I  and  II. 


TABLE  I 


i 

PROJECTILE  DATA  | 

Type 

Diameter 

Length 

Weight* 

1 

1 

\ 

1 

0.218  inch 

0.250  inch 

2.8.2  grains 

2 

0.218  inch 

O.375  inch 

27.7  grains 

1 

3 

0.299  inch 

0.341  inch 

47.3  grains 

i 

j 

i 

Material  -  434(1  steel -normalized 

1 

I 

Yield  stress  -  280,000  psi 

1 

Hardness  -  Brinell  388;  Rockwell 

CM 

1 

0 

i  8 


^approximate 


TABLE  II 


METAL  TARGET  PLATE  DATA 

Material 

Thicknesses  (inches) 

Mean  Value 

Yield  Stress  (psi) 

Hardness 

Aluminum 

7075-T6 

0.010,  0.050, 

0.090,  0.250 

74, 000 

Vickers  180 
Rockwell  A  55-5 

Titanium 

2.5  Sn  5  A1 

0.032,  0.048, 

0.076,  0.100 

117, 500 

Vickers  320 
Rockwell  A  67.75 

The  experimental  portions  of  this  study  had  two  phases. 

Phase  I  -  V..  (residual  velocity)  vs  V_  (striking  velocity)  curves  were 
■  1  X*  s 

measured  for  four  thicknesses  for  each  of  two  target  materials,  7075- T6  aluminum 
alloy  and  Ti  2.5  Sr.  5-^1  alloy,  using  a  single  type  of  projectile.  Type  1  projectile 
was  used  for  this  experimental  phase.  These  tests  provided  the  data  for  verifi¬ 
cation  and  modification  of  the  existing  perforation  models. 

Phase  II  -  Determination  of  V  vs  V  curves  for  the  thicknesses  of  one 
-  r  s 

target  material  using  two  different  types  of  projectiles.  The  Phase  II  tests 
were  conducted  using  projectile  types  2  and  3  and  0.050  and  0.090  7075-To 
aluminum  target  plates.  These  tests  provided  data  for  checking  the  modified 
perforation  models  and  for  distinguishing  between  the  perforation  models. 

The  physical  properties  of  the  target  and  projectile  materials  were 
determined.  Material  properties  determined  were  the  stress-strain  curves, 

Young's  modulus  and  hardness. 

Details  of  the  experimental  program  and  the  experimental  data  are  given 


in  Appendix  A. 


3-2  Critical  Perforation  Velocity  Measurements 


A  summary  of  the  critical  velocity  for  perforation,  V  ,  generated  on  this 
program  is  given  in  Table  III. 

TABLE  III 


TARGET  MATERIAL,  THICKNESS,  V 

’  ’  c 


Critical  Perforation  Velocity, 

r  y  (fps) 

Pro.iectile  I 

Material 

Thickness  (inches) 

Type  1 

Type  2 

Type  3 

Aluminum 

.250 

1450 

7075- T6 

.090 

750 

650 

650 

.050 

650 

550 

O 

O 

cT\ 

.010 

450  +  100 

Titanium 

.097 

1200 

5  A1  2.5  Sn 

.076 

1050 

.048 

900 

.032 

800  +  100 

*Velues  of  V  are  +  50  fps  unless  otherwise  indicated. 

0  * 


3.3  Comparison  of  Theory  and  Experiment 

The  critical  perforation  velocities  were  computed  for  three  models  of 
perforation  for  the  two  target  materials  and  the  projectile  1.  The  material 
parameters  used  in  the  computation  are  listed  in  Table  IV. 


TABLE  IV 


MATERIAL  PARAMETERS  USED  FOR  THEORETICAL  CALCULATIONS 


Projectile: 


Material 

4340  Steel 

4340  Steel 

4340  Steel 

Length 

.25" 

•  37" 

•  37" 

Radius 

.11" 

.11" 

.15" 

Density 

.283  lbs/in3 

.283  lbs/in3 

.283  lbs/in3 

aty 

180,00C  psi 

180,000  psi 

180,000  psi 

Target  Plates: 


Material 

Density 

E 

G 

O  tensile  yield  stress 

a  ,  tensile  ultimate 
tu 

h  thickness  range 


A1  7075-T6 
0.100  lbs/in3 

r 

10.3  x  10  psi 
c 

3.4  x  10  psi 
67,100  psi 
77,000  psi 

0.0"-. 25" 


Ti  5A1-2.5  Sn 

•l6l  lbs/in3 

.  6 
16  x  10  psi 

g 

6  x  10  psi 
120,000  psi 
125,000  psi 


The  predictions  for  the  critical  perforation  velocities  for  normal  impact 
of  right  circular  cylinders  of  4340  steel  into  Ti  5A1  2.5  Sn  and  7075-T6  A1  plate 
are  shown  in  Figures  4,  5,  6,  and  7>  Three  theoretical  models  of  perforation 
are  shown  on  the  graphs.  They  are: 

1.  Fracture  in  direct  wave 

2 .  Petalling  failure 

3.  Modified  shear  plug  failure 

The  predicted  critical  perforation  velocities  were  compared  with  the 
measured  critical  velocities  for  both  plate  materials  for  a  given  projectile 
(Phase  I).  The  property  of  the  plate  system  that  was  varied  was  the  plate 
thickness  (Figures  4  and  5)- 


Critical  Perforation  Velocity  (Ft/Sec) 


(Plate  Thickness )/ (Projectile  Radius) 

Figure  4  CRITICAL  PERFORATION  VELOCITY  VERSUS  PLATF  THICKNESS 
PLATE :  Ti  5A1  2 . 5  Sn 

PROJECTILE:  4340  STEEL  CYLINDER  a  =  0.109  IN.  L  =  0.25  IN 
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(Plate  Thickness )/(Projectile  Radius) 

Figure  5  CRITICAL  PERFORATION  VELOCITY  VERSUS  PLATE  THICKNESS 
PLATE:  7075 -T6  A1 

PROJECTILE:  4340  STEEL  CYLINDER  a  =  0.109  IN,  L  =  0.25  IN. 


vn 


Critical  Perforation  Velocity  (Ft/Sec) 


Figure  6  CRITICAL  PERFORATION  VELOCITY  VERSUS  PLATE  THICKNESS 
PLATE:  7075-T6  A1 

PROJECTILE:  4340  STEEL  CYLINDER  a  =  0.15  IN.  L=  0.33  IN 


0  0.5  1.0  1.5  2.0  2.5 

(Plate  Thickness  /Projectile  Radius) 


gure  7  CRITICAL  PERFORATION  VELOCITY  VERSUS  PLATE  THICKNESS 
PLATE:  7075-T6  A1 

PROJECTILE:  4340  STEEL  CYLINDER  a  =  0.109  IN.  L  =  0.375  IN 


One  of  the  theoretical  models  provided  a  good  qualitative  and  quantitative 
fit  to  the  experimentally  determined  critical  perforation  velocities  over  the 
range  of  plate  thickness  and  velocities  tested.  This  model  was  the  "modified 
shear  plug"  mode  of  failure. 

Two  of  the  theoretical  curves  used  in  combination  gave  a  reasonable 
qualitative  and  quantitative  description  of  the  critical  limit  velocity  for 
perforation.  These  modes  were  the  "petalling"  mode  for  very  thin  plates  and 
the  fracture  in  the  direct  wave  mode  for  thicker  plates,  the  critical  velocity 
being  determined  by  taking  the  lower  of  critical  velocities  predicted. 


The  Phase  II  experiments  were  determined  from  properties  of  the  modes 
of  failure.  The  "petalling"  and  "modified  shear  plug"  modes  of  failure  predict 
an  influence  of  projectile  length  (mass  per  unit  striking  area)  of  the  pro¬ 
jectile  on  the  critical  perforation  velocity.  These  curves  have  opposite 
convexities  in  the  region  of  interest.  The  critical  perforation  velocity  for 
the  combined  model  description  will  follow  either  the  fracture  in  the  direct 
wave  curve  or  the  petalling  curve.  If  the  experimental  data  is  described  by 
the  direct  wave  curve,  the  value  of  the  critical  perforation  velocity  will  be 

h  h  /  ih\2V1/2 

a  function  of  —  and  independent  of  — ,  being  proportional  to  ll  +  I— |  1 

If  the  experimental  data  is  described  by  the  petalling  mode,  the  critical 

h  i1/2/  h  iL/2  h  I  //  L  i1/2 

perforation  velocity  is  proportional  to  I—  I  (—1  =  •  If  the 

modified  shear  plug  model  is  descriptive  of  the  critical  failure  of  the  plate, 

the  behavior  of  the  critical  perforation  velocity  is  proportional  to 

/hi1/2  /ht1/2,/ hi1/2 


Near  the  intersection  of  the  direct  wave  failure  curve  and  the  petalling 
curve,  the  differences  among  these  three  theoretical  curves  is  quite  large. 
Therefore,  the  experiments  on  Phase  II  were  planned.  Additional  critical 
perforation  velocities  were  determined  for  one  of  the  plate  materials, 

7075-T6  Al.  In  one  set  of  experiments,  the  size  of  the  projectile  was  in¬ 
creased  while  retaining  geometric  similarity  to  the  smaller  projectiles 
(Projectile  3)*  The  ratio  (h/a)  was  increased  simultaneously  with  L/a. 

The  ratio  h/L  was  held  constant.  In  the  other  set  of  experiments,  the  length 
of  the  projectile  was  increased  but  the  striking  area  was  kept  the  same  as 
the  previous  projectile  (Projectile  2).  The  ratio  h/a  was  held  constant, 
while  L/a  and  consequently,  h/L  were  changed.  The  theoretical  critical 
velocities  were  determined  for  these  two  cases.  The  critical  perforation 
velocities  were  measured  for  these  cases.  Figures  6  and  7  show  the  experi¬ 
mentally  determined  values  and  the  analytical  prediction. 

The  Phase  II  experiments  indicate  that  the  "modified  shear  plug"  mode 
of  failure  gives  the  best  description  of  the  perforation  limit  velocity  for 
the  plate  materials  tested. 

One  point  remained  to  be  determined.  Was  the  effective  stress  acting 
against  the  projectile  approximated  by  the  yield  stres:  or  the  ultimate  stress? 

The  tensile  yield  stress  and  the  tensile  ultimate  stress  of  the  titanium 
and  the  aluminum  alloys  tested  are  about  the  same.  For  certain  steels  this 
is  not  the  case  and  the  experimental  determination  of  p<  ’foration  velocities 
for  steel-steel  impact  will  clarify  whether  0^  or  c s  descriptive  of  the 
impact  process.  Figure  8  shows  the  summary  of  the  NRDC Experimental  data 
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(Curtiss  1946)  on  low  velocity  steel-steel  impact  together  with  the  "modified 
shear  plug"  model.  (The  critical  velocity  for  perforation  was  modified  slightly 
for  nose  shape.  See  Vol.  II  for  details.)  The  projectiles  were  .244  cal  4340 
steel  bullets  and  the  plate  was  hardened  armor  steel  with  varying  Brincll 
hardness  (250-450). 

The  tensile  ultimate  stresses  for  armor  steel  and  naval  ship  plate  steel 
(Bridgman  1952)  Ttrera used  to  compute  the  critical  velocity  curves.  The  modified 
shear  plug  model  was  used.  Good  agreement  using  tensile  ultimate  stress  on 
armor  steel  was  obtained  for  plate  thickness-projectile  radius  ratios  from 


0  to  12. 


SECTION  4 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  main  goal  and  objective  of  this  research  project  was  the  determination 
of  the  descriptive  model  of  failure  of  thin  metallic  plates  under  ballistic  im¬ 
pact  by  metallic  projectiles.  The  experimental  discrimination  and  verification 
of  the  model  of  failure  was  accomplished. 

A  review  of  previous  analytical  research  on  the  determination  of  the  stresses 
in  the  plate  due  to  impact  led  to  the  delineation  of  several  models  of  failure  of 
the  plates.  For  each  model  of  failure, analytic  expressions  for  the  critical  per¬ 
foration  velocity  and  for  the  residual  velocity  versus  striking  velocity  were 
formulated  as  a  function  of  plate  thickness,  projectile  geometry,  and  plate  and 
projectile  properties.  Calculations  of  ballistic  parameters,  in  particular  the 
critical  perforation  velocity,  were  made  for  the  ballistic  impact  of  two  thin 
metallic  plates  composed  of  an  aluminum  alloy  and  a  titanium  alloy,  respec¬ 
tively,  by  right  circular  cylinders  of  steel. 

The  critical  perforation  velocity  and  the  residual  velocity  versus  striking 
velocity  curves  were  experimentally  determined  for  the  normal  impact  of  right 
circular  cylinders  of  steel  into  these  two  alloy  plates.  The  experimentally 
determined  critical  perforation  velocities  and  residual  velocities  agreed 
very  well  with  one  of  the  theoretical  models  for  failure  of  the  plate  over  the 
range  of  plate  thicknesses  and  impact  velocities  tested.  The  plate  thickness 
varied  up  to  0.25"  and  the  normal  impact  velocity  range  covered  velocities  to 
4000  ft/sec. 
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The  particular  model  of  plate  failure  that  best  describes  the  ballistic 
performance  of  a  thin  metallic  plate  under  the  ballistic  attack  of  cylindrical 
projectiles  is  the  "modified  shear  plug"  model  of  failure.  In  this  model  of 
plate  failure,  the  impacting  projectile  pushes  the  material  in  the  plate  out 
of  the  way.  The  resisting  force  of  the  plate  against  the  projectile  motion 
is  composed  of  the  inertial  resistance  of  the  plate  material  and  the  compressive 
ultimate  stress  of  the  plate  in  unconfined  compression.  The  critical  perforation 
velocity  and  tne  residual  velocities  are  expressible  in  closed  form  for  the 
cylindrical  projectile.  The  computations  are  easily  extended  to  the  cases  of 
normal  impact  of  any  axisymmetric  projectile.  Representative  calculations  for 
projectiles  in  this  class  are  given  in  the  second  volume  of  this  report. 

The  critical  velocity  for  perforation  for  a  cylindrical  projectile  is  a 
function  of  the  compressive  ultimate  stress  of  the  plate  material,  the  density 
of  the  plate  material,  the  density  of  the  projectile  and  the  plate-  thickness 
to-projectile  length  ratio.  The  critical  velocity  for  perforation  is  propor¬ 
tional  to  the  ultimate  stress  of  the  plate  and  inversely  proportional  to  the 
plate  density.  The  critical  velocity  is  a  monotonically  increasing  function 
(almost  exponential)  of  the  ratio  of  the  mass-per-unit  area  of  the  plate  to 
mass«.penrunit  area  of  the  projectile. 

After  the  critical  perforation  velocity  has  been  surpassed,  the  residual 
velocity  to  striking  velocity  ratio  rapidly  increases  with  increasing  value 
of  the  striking  velocity.  This  ratio  soon  reaches  a  limiting  value  less  than 
unity. 
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One  of  the  ultimate  goals  of  this  research  program  was  to  give  information 
which  will  increase  the  ballistic  performance  of  light  armor  systems.  The 
second  volume  of  this  report  applies  two  of  these  models  of  perforation  to 
evaluate  ballistic  performance  of  light  metallic  plates  by  cylinders  and 
military  projectiles.  Recommendations  on  the  improvement  of  armor  systems 
is  given  in  that  volume. 


APPENDIX  A 

EXPERIMENTAL  DETERMINATION  OF  THE 
RESIDUAL  VELOCITIES  FOR  NORMAL  IMPACT  OF 
RIGHT  CIRCULAR  CYLINDERS  INTO  THIN  METALLIC  PLATES 
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A.l  Introduction 


Experimental  determination  of  the  ballistic  performance  of  thin  metallic 
alloy  plates  subject  to  normal  impact  by  right  circular  steel  cylinders  was 
made.  The  details  of  these  experiments  are  given  in  this  appendix. 

The  experiments  determined  the  residual  velocities  and  the  critical  limit 
perforation  velocities  for  the  normal  impact  of  4340  steel  cylinders  into 
selected  titanium  and  aluminum  alloys.  This  experimental  program  coincided 
with  the  readily  computable  cases  generated  by  the  theoretical  and  analytical 
phases  of  the  perforation  problem.  The  experiments  were  designed  to  check 
the  velocity  of  the  various  models  of  perforation  and  to  discriminate  among 
the  various  models. 

A. 2  Experimental  Setup 

A. 2.1  Apparatus 

The  experimental  setup  consisted  of  a  gun  fixture,  a  test  fixture  which 
held  the  target  plate,  a  ballistic  pendulum  and  several  velocity  screen  mounts. 
It  is  shown  schematically  in  Figure  A-l. 

The  gun  fixture  and  the  test  fixture  are  separated  by  about  6  feet. 

Along  the  projectile  flight  path,  a  series  of  velocity  screens  (paper  sheets 
with  a  continuous,  printed  silver  circuit)  and  a  ballistic  pendulum  are  placed. 
Four  screens  are  placed  between  the  gun  muzzle  and  target  plate.  Two  screens 
are  placed  between  the  target  and  the  ballistic  pendulum, which  is  immediately 
behind  the  target  fixture.  Three  average  velocities  are  obtained  in  front  of 
the  target  plate  and  1  velocity  is  obtained  behind  for  the  screen  measurements. 


Plat 


Figure  A-l  DIAGRAM  OF  EXPERIMENTAL  SETUP 


The  ballistic  pendulum  is  about  three  meters  in  suspended  length  and  its  mass 
may  be  varied  fran  about  two  kilograms  to  nine  kilograms. 

Target  plates  varying  in  width  from  about  one  foot  to  three  feet  are 
easily  acccmnodated  in  the  plate  holder.  Two  types  of  projectile  launchers 
were  used:  a  caliber  .220  Swift  type  for  projectiles  from  0.2180  to  0.2185 
inches  in  diameter,  and  a  caliber  .350  Magnum  type  for  projectiles  from  0.2990 
to  0.2995  inches  in  diameter.  Both  of  the  launchers  have  smooth  bore  barrels. 

A. 2. 2  Experimental  Procedure 

Measurements  of  the  velocity  of  the  projectile  before  and  after  it  per¬ 
forated  the  target  plate  were  made  for  a  sequence  of  initial  velocities. 

The  sequence  of  the  experimental  shots  was  the  same  for  all  target  plates. 

The  sequence  consisted  of  eight  pairs  of  shots  in  three  sets.  The  first  set  of 

three  pairs  of  shots  was  fired  at  a  velocity  somewhat  above  the  knee  of  the 

analytically  predicted  V  versus  V  curve  (See  Figure  A-2).  The  second  set  con- 

r  s 

sisted  of  one  pair  at  the  high  end  of  the  curve  and  one  pair  near  the  knee  of 
the  curve.  The  third  set  of  three  pairs  were  near  the  value  of  the  curve. 

The  velocities  used  in  sets  2  and  3  were  dependent  upon  the  previously  obtained 
results.  Thus,  for  each  plate  tested,  a  versus  Vg  curve  was  obtained  for 
comparison  with  the  analytically  determined  curves. 

Visual  inspection  and  metallurgical  examination  of  the  plate  in  the 
neighborhood  of  the  perforation  were  made  as  necessary  to  assist  in  the  deter¬ 
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mination  of  the  mode  of  failure. 


Residual  Velocity 


Striking  Velocity 


Figure  A-2  GROUPING  OF  TEST  SHOTS  FOR  DETERMINING  RESIDUAL  VELOCITY 
VERSUS  STRIKING  VELOCITY  CURVE 


Several  mechanical  properties  of  the  target  plate  were  also  determined  as 
a  part  of  the  experimental  portion  of  this  study. 

A. 2. 3  Special  Problems  and  Solutions 

During  the  initial  stages  of  the  experimental  portion  of  this  study,  several 
problems  were  encountered. 

One  problem  encountered  envolved  the  measurement  of  the  residual  velocity 
with  a  decade  counter.  During  the  early  phase  of  this  study,  a  determination  of 
the  residual  velocity  by  means  of  a  pair  of  velocity  screens  spaced  five  centi¬ 
meters  apart,  close  behind  the  target  plate,  was  attempted.  These  screens  were 
connected  to  a  sensitive  decade  counter. 

On  the  initial  tests,  the  residual  velocity  measurements  obtained  from  the 
counter  differed  significantly  from  the  ballistic  pendulum  measurements,  partic¬ 
ularly  at  higher  velocities  (see  decade  counter  measurements  on  Figure  A-3) .  This 
discrepancy  was  thought  due  to  small  particles  coming  off  the  plate  prior  to  the 
plug  and/or  the  projectile  coming  through  the  plate.  Th<|se  small  particles 
were  traveling  at  velocities  near  the  striking  velocity,  but  had  energy  only 
sufficient  to  interrupt  the  first  velocity  screen,  starting  the  decade  counter. 

The  counter  was  stopped  when  a  more  massive  particle,  such  as  the  projectile 
or  plug,  interrupted  the  second  screen.  Two  attempts  were  made  to  correct  this 
situation.  First,  a  spall  shield  was  used  in  front  of  the  first  velocity  screen 
after  the  target  plate.  This  was  a  wooden  plate  with  a  hole  in  it  which  would 
permit  the  passage  of  only  those  particles  leaving  the  plate  almost  normally. 

This  did  not  solve  the  problem.  Figure  A-5  shows  the  residual  velocity  measure¬ 
ments  with  the  spall  shield  in  place.  Second,  the  velocity  screens  were  moved 


away  from  the  plate.  The  "punchout"  material  was  forced  to  pass  through  several 
sheets  of  paper  before  encountering  the  velocity  screens.  This  also  failed  to 
solve  the  problem.  Figure  A-’+  shows  the  velocity  measurements  with  the  paper 
sheets  in  place.  Interposing  several  thicknesses  of  paper  was  used  while  firing 
at  the  0.090"  aluminum  plate.  The  intent  was  to  stop  all  particles  other  than 
the  projectile  and  plug.  The  effect  of  the  paper  was  to  create  a  greater  di¬ 
vergence  between  the  measured  ballistic-pendulum  and  decade-counter  velocities. 

This  result  indicated  that  the  particles  had  not  been  stopped,  even  when  an 
excessive  number  of  paper  sheets  were  used  (l6  sheets).  While  the  use  of  such 
a  quantity  of  paper  might  not  affect  the  ballistic  pendulum  reading  appreciably 
at  higher  residual  velocities,  it  most  certainly  would  affect  the  readings  at 
the  lower  residual  velocities. 

It  appeared  unwarranted  to  pursue  this  problem  further,  and  this  velocity 
measurement  system  was  discontinued.  Another  factor  leading  to  the  discontinuation 
was  that  about  10  percent  of  the  counter  readings  were  obviously  spurious.  The 
cause  of  the  obviously  spurious  results  was  not  found,  but  itv6s  possibly  due 
to  the  high  sensitivity  of  the  counter  and  the  proximity  of  the  other  electronic 
equipment . 

Another  problem  expected  and  encountered,  was  the  determination  of  the  mass 
of  the  material  affecting  the  ballistic  pendulum.  In  the  initial  series  of  tests, 
the  mass  affecting  the  pendulum  could  not  be  determined  with  as  high  a  degree  of 
accuracy  as  had  been  anticipated.  Therefore,  the  ballistic  pendulum  was  modified 
so  that  a  large  percentage  of  the  mass  affecting  the  pendulum  could  be  recovered 


and  its  weight  determined. 


Another  minor  problem  encountered  was  that  of  the  projectile  tumbling  in 
flight  before  striking  the  target  plate.  There***-  a  high  tendency  for  the  short 
cylindrical  projectile  to  tumble  at  muzzle  velocities  above  36OO  fps,  and  a 
slight  tendency  for  tumbling  to  occur  in  the  500  fps  to  1000  fps  muzzle  velocity 
range.  The  tumbling  of  the  projectile  was  apparently  confined  to  these  areas 
and  was  relatively  intermittent  even  within  these  areas.  No  special  attempt  was 
made  to  correct  the  situation. 

On  the  initial  set  of  shots,  some  of  the  projectiles  did  not  strike  the 
target  plate  normally.  A  few  projectiles  were  canted  as  much  as  15  degrees 
when  striking  the  target  plate.  The  misalignment  of  this  order,  or  less,  does 
not  change  the  qualitative  features  of  the  deformation  and  perforation  processes. 
Its  quantitative  effect  in  this  velocity  regime  is  assumed  to  be  small. 

A. 2. 4  Accuracy  of  Results 

The  experimental  results  are  sufficiently  accurate  for  the  purpose  of  this 
study.  A  maximum  error  of  +  10  percent  in  the  residual  velocities  could  exist 
at  low  residual  velocity  (<  500  fps),  but  this  error  decreases  as  the  residual 
velocity  increases  (at  of  about  2000  fps,  the  error  is  +  2  percent).  The 
maximum  absolute  error  in  the  residual  velocity  expected  in  our  experiments  is 
on  the  order  of  +  30  feet  per  second. 

The  accuracy  of  the  test  results  is  primarily  fixed  by  the  accuracy  of  the 
determination  of  the  mass  of  the  pendulum  and  the  measurement  of  the  pendulum 
swing  .stance. 


The  distance  through  which  the  pendulum  moves  is  easily  determined  within 
0.05  centimeters.  For  the  majority  of  the  tests  which  resulted  in  residual 
velocities  of  less  than  500  fps,  the  swing  distance  of  the  pendulum  was  more 
than  1  centimeter.  The  maximum  error  of  the  pendulum  swing  measurement  was 
+  5  percent.  The  range  of  the  residual  velocity  was  known  and  a  relatively 
large  movement  of  the  pendulum  (hence  a  smaller  error)  was  obtained  by  proper 
variation  of  the  pendulum's  mass. 

The  determination  of  the  mass  affecting  the  pendulum  is  more  difficult 
since  it  is  not  practical  to  recover  every  bit  of  the  mass  from  the  pendulum. 
Fortunately,  this  is  not  required.  If  all  mass  other  than  the  projectile  were 
disregarded,  the  error  in  the  residual  velocity  would  be  less  then  +  10  percent 
in  all  the  plates  tested  with  the  exception  of  tte  0.25-inch  aluminum  and  0.1-inch 
titanium  plates.  For  these  two  plates,  the  error  incurred  by  disregarding  the 
mass  other  than  the  projectile  could  run  as  high  as  +  25  percent.  The  recovery 
of  any  relatively  large  mass  other  than  the  projectile,  such  as  a  plug,  con¬ 
siderably  reduced  these  error  figures. 

The  errors  in  the  determination  of  the  residual  velocity  were  less  than 
+  5  percent  for  V  less  than  1000  fps,  and  less  than  +  3  percent  for  V  greater 

—  r  a. 

than  1000  fps, 

A. 2. 5  Reduction  of  the  Data  and  the  Determination  of  the  Critical 
Perforation  Velocity 

The  experimental  data  furnished  are  a  sequence  of  striking  velocities  (V  ) 

s 

and  the  corresponding  residual  velocity  after  perforation  (Vr).  If  the 


projectile  did  not  perforate  the  plate,  the  residual  velocity  was  zero.  One  of 

the  main  problems  of  the  experimental  procedure  is  the  determination  of  the 

critical  perforation  velocity,  V  .  This  velocity  is  the  striking  velocity 

c 

below  which  there  is  no  perforation  and  above  which  there  is  perforation. 

The  critical  perforation  velocity  cannot  be  precisely  measured.  A  number 
of  factors  with  random  variations,  such  as  canting  or  nonuniformity  of  plate 
strength,  limit  the  sharpness  of  the  transition  from  nonperforating  impact  to 
perforating  impact.  A  preliminary  estimate  of  the  critical  perforation  velocity 
tiras  obtained  by  taking  the  highest  experimental  velocity  where  no  perforation 
occurred.  This  is  a  lower  bound  to  the  critical  perforation  velocity.  The  dif¬ 
ficulty  in  obtaining  precisely  the  striking  velocity  near  the  critical  velocity 
requires  an  extrapolation  from  the  residual  velocities  at  higher  striking 
velocities . 

The  analytical  portion  of  the  study  indicates  that  the  residual  velocity- 
striking  velocity  curve  is  of  the  form  for  a  given  plate  and  projectile. 


where  and  are  constants . 

A'  least  square  fit  of  the  residual  velocity  and  the  striking  velocity 
obtained  the  constants  and  C^.  Tne  critical  perforation  velocity,  Vq, 


is  the  striking  velocity  at  which  the  residual  velocity  just  goes  to  zero.  Thus, 


V  = 
c 

The  computation  of  the  critical  perforation  velocity  by  the  direct  deter¬ 
mination  of  the  constants  was  very  inaccurate.  The  following  procedure  was 
applied  to  determine  V  .  The  values  of  V^/V^  were  plotted  versus  V^.  The 

limiting  value  of  V  /V  as  V  becomes  large  was  easily  determined.  Denote  this 

r  s  s 

limiting  value  by  (V  /V  )  ».  The  striking  velocity  at  which  the  ratio  V  /V 

c'  s  r  s 

attains  one-half  of  the  limiting  ratio  was  also  determined.  Denote  this  velocity 
by  V^.  Using  the  analytical  form  of  the  residual  velocity  curve, 

V  =  0.89  V. 

c  1 

A. 3  Experimental  Results 

A. 3.1  Phasing  of  Test  and  General  Data 

The  experiments  obtained  the  residual  velocity  vs.  striking  velocity  curves 
for  use  in  checking  and  modifying  various  penetration  models.  Three  types  of 
projectiles  were  shot  at  various  thicknesses  of  the  types  of  metal  target 
plates. 

Properties  of  the  three  types  of  projectiles  and  the  metal  plates  used 
in  the  experimental  phases  of  this  study  are  given  in  Tables  I  and  II. 


TABLE  I 


PROJECTILE  DATA 

type 

Diameter 

Length 

Weight* 

1 

0.218  inch 

0.250  inch 

18.3  grains 

2 

0.218  inch 

0.375  inch 

27.7  grains 

3 

0.299  inch 

0.34l  inch 

47.3  grains 

Material  -  4340  steel-normalized 

Yield  Stress  -  140,000  psi 

Hardness  -  Brinell  388  Rockwell 

C-42 

^approximate 


TABLE  II 


METAL  TARGET  PLATE  DATA 

Material 

Thickness  (inches) 

Yield  Stress  (psi) 

Hardness 

Aluminum 

0.010,  0.050, 

74,000 

Vickers  180 

0.090,  0.250 

Rockwell  A  55-5 

Titanium 

0.032,  0.048, 

117, 500 

Vickers  320 

0.076,  0.097 

Rockwell  A  67.75 

The  experimental  portions  of  this  study  were  conducted  in  two  phases: 

Phase  I  -  Residual  velocity  (V_ )  versus  striking  velocity  (V_)  curves 
were  determined  for  four  thicknesses  for  each  of  two  target  materials  using 
a  single  type  of  projectile.  These  tests  provided  data  for  verification  and 
modification  of  the  existing  penetration  models. 


Phase  II  -  V  vs  curves  were  determined,  for  two  thicknesses  of  one  target 

■  ■■■■--■"  r  s  ° 

material  using  two  different  types  of  projectiles.  These  tests  provided  data  for 
checking  the  modified  penetration  models. 

An  additional  phase  of  the  experimental  study  was  the  determination  of  the 
physical  properties  of  the  target  and  projectile  materials.  Material  pr-  'rties 
obtained  were  stress-strain  curves.  Young's  modulus,  and  the  tensile  yield  stress. 

A. 3. 2  Critical  Velocity  and  Residual  Velocity  Measurements 

A  summary  of  the  critical  perforation  velocities,  V  ,  generated  on  this 

c 

program  is  given  in  Table  III. 

Phase  I  Tests  -  V  vs  V  curves  were  determined  for  four  thicknesses  of 

r  s 

7075-T6  aluminum  target  material  and  for  four  thicknesses  of  titanium  2.5  Sn 
5  A1  aluminum  alloy  target  material.  Type  1  projectiles  were  used  for  this 
experimental  phase. 


TABLE  III 


TARGET  MATERIAL 

,  THICKNESS, 

V 

c 

Critical  Perforation  Velocity 

,  Vc*  (fps) 

Projectile 

Material 

Thickness  (inches) 

1 

2 

3 

Aluminum 

.250 

1450 

7075-T6 

.090 

750 

650 

650 

.050 

650 

550 

500 

.010 

500  +  100 

Titanium 

.097 

1200 

5  A1  2.5  Gn 

076 

1050 

.048 

900 

.032 

800  +  100 

*Valuos  of  V  are  +  50  fps  unless  otherwise  indicated. 


'•***■**— r.v  ^  '»i«i  *  *"■ 
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Figures  3  through  10  present  the  data  obtained  in  Phase  I  and  show  both  the 

experimental  points  and  an  approximate  V  vs  V  curve.  Figures  11  and  12  show 

r  s 

the  Vr  vs  Vs  curves  of  the  two  tested  target  materials  as  they  appear  as  a 
family. 

Phase  II  Tests  -  Phase  II  tests  were  conducted  using  projectile  types  2  and 
3  and  0.050  and  0.090  7075-T6  aluminum  target  plates. 

Figures  13,  14,  15,  and  16  present  the  data  obtained,  showing  both  the 

experimental  points  and  an  approximate  V  vs  V  curve.  The  family  of  V  vs  V 

t  s  r  s 

curves  of  the  three  types  of  projectiles  for  each  of  the  two  plate  thicknesses 
tested  are  shown  in  Figure  17  and  18. 

A  summary  of  the  material  properties  tests  -ia  given  in  Table  IV  and  a 
typical  set  of  data  is  shown  plotted  in  Figure  19, 
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Striking  Velocity  -  feet  per  second 

Figure  A-3  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  FOR  0.250"  7075-T6  ALUMINUM  TARGET  MATERIAL 


Projectile  -  4340  Steel  -  Normalized 
Diameter  -  0.218" 


Striking  Velocity  -  feet  per  second 

Figure  A-4  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  FOR  0.090"  7075-T6  ALUMINUM  TARGET  MATERIAL 


Figure  A- 5  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  FOR  0.050"  7075-T6 


Projectile  -  4340  Steel  -  Normalized 

Diameter  -  0.218" 

Length  -  0.250" 

Weight  -  l8.3  Grains 

A  -  Experimental  Pointr, 
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Figure  A- 6  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  FOR  0.010”  7075-T6  ALUMINUM  TARGET  MATERIAL 
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Striking  Velocity  -  feet  per  second 

Figure  A- 12  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  AS  A  FUNCTION  OF  TARGET  THICKNESS 
FOR  TITANIUM  ALLOY  TARGET  METER IAL 


Projectile  -  4340  Steel  -  Normalized 

Diameter  -  0.218” 

Length  -  0.375" 

Weight  -  27.6  Grains 
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Figure  A- 13  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  FOR  0.050"  7075-T6  ALUMINUM  TARGET  MATERIAL 
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Figure  A-l4  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  FOR  0.090"  7075-T6  ALUMINUM  TARGET  MATERIAL 


Projectile  -  4340  Steel  -  Normalized 

Diameter  -  0.299” 

Length  -  0.341" 

Weight  -  47.3  Grains 

A  -  Experimental  Point 
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Figure  A- 15  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  FOR  0.050"  7075-T6  ALUMINUM  TARGET  MATERIAL 
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Figure  A- 17  RESIDUAL  VELOCITY  VS.  STRIKING  VELOCITY  AS  A  FUNCTION  CF  PROJECTILE  GEOMETRY  FOR 
0.050  INCH  THICK  7075-T6  ALUMINUM 


TABLE  IV 


MATERIA,  PROPERTIES 

Material:  7075- T6  Aluminum 

Yeung's 

Yield 

Hardness 

Ultimate 

Thickness 

Modulus 

Stress 

Rockwell 

Stress* 

(inches) 

(psi) 

(psi) 

Vickers 

A-Scale 

(psi) 

0.250 

10.2  x  10 

75,500 

190 

56.0 

83, 500 

0.090 

10.1  x  106 

74, 500 

195 

57.0 

84,500 

0.050 

10.0  x  106 

74,000 

188 

55.5 

83,000 

0.010 

9.8  x  106 

72,000 

177 

(53.5)** 

78,000 

Mean  Values 

74, 000 

188 

55.5 

82,000 

j  Material:  Titanium  Alloy 

Ti  5  A1  2.5 

Sn 

0.100 

16.5  X  106 

116, 500 

(317) 

66.5 

120,000 

0.070 

£ 

16.0  x  10 

115, 500 

(325) 

67.0 

122,000 

0.049 

16.8  x  106 

122,000 

(309) 

66.0 

130,000 

0.033 

16.1  x  106 

118,000  , 

(334) 

67.5 

132,500 

Mean  Values 

117, 500 

(320) 

66.75 

126, 500 

*  Estimate  based  upon  final  tensile  test  load. 

**  Numbers  in  ()  represent  equivalents,  not  experimental  data. 
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APPENDIX  B 

PHYSICS  OF  THE  PERFORATION  OF  THIN 
METALLIC  PLATES  BY  PROJECTILE  IMPACT 
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B.l  Ma.jor  Problem  Areas 

This  appendix  presents  a  summary  of  the  analytical  and  the  theoretical 
studies  on  penetration  and  perforation  of  thin  metallic  plates  by  normal  im¬ 
pact  of  projectiles.  For  the  most  part,  this  appendix  is  a  review  of  the 
theoretical  models,  assumptions,  and  results  of  the  previous  reports  (Fugelso, 
Arentz,  and  Foczatek  1961,  and  Fugelso,  Arentz,  and  Davidson  1962).  Some  cor¬ 
rections  have  been  made  where  necessary  and  sane  points  have  been  more  fully 
described.  For  the  mathematical  detail  the  reader  is  referred  to  the  previous 
studies. 

Consider  a  right  circular  cylindrical  projectile  of  radius,  a,  length, 

L,  and  mass,  M,  approaching  a  plate  of  uniform  thickness,  h.  The  projectile 

is  oriented  so  that  its  axis  is  perpendicular  to  the  plate.  The  projectile 

approaches  the  plate  with  a  velocity,  V  ,  which  is  perpendicular  to  the  surface 

s 

of  the  plate.  There  is  no  lateral  velocity  component.  When  the  projectile 
strikes  the  plate,  the  entire  flat  end  of  the  cylinder  will  strike  the  flat 
surface  of  the  plate.  The  density  of  the  plate  is  p ^  and  the  density  of  the 
projectile  is  PQ  (See  Figure  l).  The  question  that  must  be  answered  is, 

"Will  the  projectile  perforate  the  plate?" 

Hopkins  and  Kolsky  (i960)  presented  a  decomposition  of  impact  phenomena 
based  on  the  predominant  mechanism  of  deformation.  Regions  wherein  the  de¬ 
formation  was  predominantly  elastic  deformation,  plastic  deformation  or 
hydrodynamic  flow  were  estimated  on  the  basis  of  the  stress  generated  at 
the  surface.  The  impact  stresses  are  expressed  in  terms  of  the  impact 
velocity  for  one-dimensional  impact. 
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Figure  B-l  GEOMETRY  OF  IMPACT  AND  DEFINITION  OF  TERMS 


{ 


The  upper  hound  on  the  impact  velocity  for  the  elastic  region,  which  in 
turn  is  the  lower  bound  on  the  impact  velocity  for  the  plastic  region,  is  the 
impact  velocity  which  generates  the  static  yield  stress  in  the  plate. 


O 

v..  =  Elastic-Plastic  Transition  Velocity 

■L  p  C 

The  upper  "bound  on  the  impact  velocity  for  the  region  in  which  plastic 
deformation  is  dominant  is  the  velocity  at  which  the  hydrodynamic  decelerating 
force  on  the  projectile  is  greater  than  the  yield  strength  of  the  plate. 


Plastic-Hydrodynamic  Transition  Velocity 


Typical  values  for  the  transition  velocities  for  one-dimensional  (uni¬ 
axial)  impact  are  given  in  Table  B-l.  The  effects  of  each  mode  of  deformation 
persist  throughout  the  entire  range  of  impact  velocities.  Thus  a  more  detailed 
study  of  penetration  was  warranted  than  was  carried  out  in  the  previous  study. 


TABLE  B-l 


TRANSITION  VELOCITIES  FOR  IMPACT  (FT/SEC)* 


Transition 

Velocities 

Steel 

Duralum 

A1 

Cu 

Pb 

V1 

150 

200 

^3 

CVJ 

7 

V2 

1180 

1410 

o\ 

ro 

v/i 

425 

130 

*Projectile  and  plate  materials  are  the  same. 


To  determine  the  plate's  failure  under  ballistic  impact,  several  problems 
had  to  be  solved.  The  dynamic  problem  of  the  stress  wave  propagation  in  the 
plate  had  to  be  solved.  After  the  stresses  had  been  determined,  the  areas  of 


the  plate  where  failure  was  likely  to  occur  could  be  delineated.  Specifically, 
the  following  stress  wave  propagation  problems  were  solved: 

1.  The  transient  stress  wave  propagation  in  a  plate  system  under 
the  normal  impact  of  a  cylindrical  projectile  was  solved  analytically 
and  numerically  for: 

a.  Single  layer  of  an  isotropic  metal 

b.  Double  layer  of  an  isotropic  metal 

2.  An  approximate  solution  for  the  propagation  of  large  trans¬ 
verse  deflection  of  the  plate  was  established  numerically,  including  the 
effects  of  plastic  flow. 

It  was  found  that  several  properties  of  material  deformation  under  rapid 
transient  loading  had  to  be  defined.  Therefore,  a  phase  of  the  prior  effort 
was  devoted  to  material  properties  at  high  strain  rates.  Specifically,  the 
following  problem  areas  were  investigated: 

1.  A  theory  for  macroscopic  plastic  flow  under  transient  leading 
conditions  was  developed  using  the  properties  of  inter- and  intra- crystalline 
imperfections  (Fugelso  1962).  These  results  were  correlated  with  the  experi¬ 
mental  data  available  on  the  high  strain. 

2.  A  theory  for  fracture  of  metals  under  transient  loading  was 
developed,  utilizing  elasticity  and  the  flow  theory  developed  above, 
and  compared  with  existing  data  on  fracture  of  metals  at  high  stress 


The  solutions  to  the  stress  wave  propagation  problems  delineated  several 
regions  of  the  plate  where  the  stresses  might  cause  fracture.  Each  possible 
mode  of  fracture  of  the  plate  in  response  to  the  stresses  defined  a  model  of 
failure  of  the  plate.  Each  model  of  failure  (l)  required  certain  levels  of 
impact  stresses  to  activate  the  mode  of  fracture  and  (2)  generated  a  re¬ 
tarding  force  on  the  projectile. 

After  the  stress  wave  propagation  was  solved  and  the  models  of  failure 
delineated,  the  deceleration  of  the  projectile  was  computed,  using  the  appro¬ 
priate  interface  stresses  to  oppose  and  resist  the  projectile’s  motion. 

From  these  solutions,  predictions  for  the  critical  perforation  velocities 
and  the  residual  velocity  versus  striking  velocity  were  derived. 

The  next  section  reviews  the  problem  of  the  determination  of  the  stress 
wave  propagation,  indicates  the  specific  problems  solved,  lists  the 
assumptions  involved,  and  presents  the  solutions  and  the  results  of  the 
calculations. 

Section  B-3  defines  the  modes  of  failure  of  the  plate  and  summarizes  the 
derivations  of  the  ballistic  performance  parameters. 

B.2  Stress  Wave  Propagation  in  the  Plate 

B.2.1  Introduction 

The  propagation  of  stresses  into  the  plate  due  to  the  impact  stresses 
generated  at  the  contact  interface  were  solved  by  two  approximate  methods. 
Each  approximation  was  designed  to  demonstrate  certain  qualitative  and 
quantitative  features  of  the  transient  stress  distribution  in  the  plate. 


The  first  problem  solved  was  the  propagation  of  a  stress  wave  propagating  into 
a  plate  which  was  a  linear  elastic  medium.  The  second  problem  was  the  propa¬ 
gation  of  transverse  deflections  of  a  thin  plate  composed  of  material  whose 
description  included  plastic  deformation. 

B.2.2  Linear  Elastic  Wave  Propagation 

The  first  group  of  problems  solved  was  the  propagation  of  a  linear  elastic 
stress  wave.  These  problems  were  solved  to  determine  the  qualitative  description 
of  the  various  stress  components  throughout  the  plate  as  a  function  of  position 
and  time,  and  to  obtain  a  quantitative  estimate  of  these  stresses. 


The  equations  of  motion  for  axisymmetric  deformation  are: 


+ 


Vae 


where  a  ,  oa,  a  >  T  are  the  stresses  in  the  plate 
x  u  z  rz 

u^,  uz  are  the  radial  and  vertical  displacement 
p  is  the  density  of  the  plate 
r,  z  are  the  radial  and  vertical  coordinates 
t  is  the  time 


71 


The  linear  elastic  stress-strain  relationships  used  in  the  computation 


are 


a 

z 


bn  u 

(x  +2 m)  ~  +  + 

or 

du  u 

( X  +  2jJ)  — -  +  X(— £  + 
dz 


u  du 

(x  +  2/i)  -r  +  x{— ^  + 

dr 

du  du 

m(—  +  — ) 

dz  dr 


dz 

dr 

K 

dz 


) 

) 

) 


where  X,j u  are  Lame's  constants 

The  boundary  condition  at  the  bullet-plate  interface  was  chosen  to  be  a 
uniform  vertical  stress  acting  over  the  contact  area.  While  this  is  not  the 
exact  distribution  of  stress  expected  there,  it  is  close  enough  to  determine 
the  physics  involved  in  the  process  of  deformation  and  to  approximate  the 
perforation  parameters.  The  other  exterior  surfaces  are  free  from  all 
tractions . 

Formal  solutions  for  the  equations  of  motion  were  obtained  by  application 
of  Fourier-Hankel  integral  transforms.  These  equations  were  solved  by 
partially  analytical  and  partially  numerical  techniques  for  the  following 
three  problems: 

1.  Elastic  Stress  Waves  in  a  Half  Space:  The  stress  wave  propagation 
in  an  elastic  half  space  under  the  action  of  a  vertical  stress  applied  uni¬ 
formly  over  a  circular  region  was  solved.  Numerical  solutions  for  the 
stresses,  displacements  and  velocities  throughout  the  region  were  computed. 


Figures  B-2  and  B-3  show  the  progression  of  the  vertical  stress  into  the 
half  space  at  various  times.  Figure  B-4  shows  a  profile  of  the  vertical  stress 
on  the  axis  of  symmetry.  Figure  B-5  shows  the  vertical  and  radial  stresses  on 
the  axis  of  symmetry. 

Quantitative  and  qualitative  features  of  interest  are: 

a.  The  stress  jump  at  the  wave  front  under  the  bullet  is  almost 
uniform  with  r.  The  magnitude  of  the  jump  at  the  wave  front  decays  as 

a  =  a  °  ( l  +  (-)2)"1/2 
z  z  '  va  ' 

where  0 z°  is  the  uniform  stress  acting  on  the  surface 

b.  This  stress  spike  near  the  wave  front  decays  very  rapidly 
behind  the  shock  front. 

c.  The  given  static  solution  is  established  (to  within  5$)  after 
the  passage  of  the  shear  wave  front.  The  magnitude  of  the  static  part  here  is 

°,s  -  a  *  (ffr1 

2.  Elastic  Stress  Waves  in  Plates:  The  vertical  and  radial  stresses 
on  the  axis  of  symmetry  for  an  elastic  plate  ir.  response  to  suddenly  applied, 
uniformly  distributed  vertical  pressure  were  calculated.  Typical  values  for 
the  vertical  and  radial  stress  for  a  moderately  thick  plate  (—  =  0.8)  for  two 

cl 

different  plate  materials  are  shown  in  Figures  B-6  and  B-7. 

A  large  compressive  vertical  stress  is  generated  in  the  spike  as  the 
wave  travels  across  the  plate.  The  vertical  stress  remains  compressive  in 
the  first  wave  after  the  spike  has  decayed.  The  vertical  stress  is  reflected 
with  a  reversed  sign  at  the  frej  surface.  The  combined  vertical  stress  becomes 
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z/a 


Figure  B-2  CONTOURS  OF  VERTICAL  STRESS  (7T  a  )  UNDER  THE 

z 

BULLET  AT  VARIOUS  TIMES  AFTER  IMPACT  (POSITIVE 
SIGN  IS  COMPRESSIVE) 


Figure  B-4  VERTICAL  STRESS  VERSUS  DEPTH 
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Figure  B»?  RADIAL  A?D)  VERTICAL  STRESS  ON  THE  AXIS  OF  SYMMETRY 


B -6 


tensile  a  little  way  from  the  hack  or  free  surface.  It  is  at  this  point  where 
failure  by  spall  will  occur.  Failure  of  the  plate  material  could  also  occur  in 
the  first  passage  of  the  stress  wave  by  the  fact  that  the  compressive  stresses 
are  large. 

The  radial  stress  in  the  spike  is  compressive.  After  the  spike  has  passed, 
the  radial  stress  in  the  front  wave  becomes  tensile.  The  radial  stress  in  the 
reflected  stress  spike  is  tensile.  The  combined  radial  stress  system  adds  two 
tensile  stresses. 

After  the  complete  passage  of  the  reflected  wave  has  passed  and,  for  sub¬ 
sequent  ref lections, the  vertical  stress  in  the  plate  is  quite  small. 

3.  Elastic  Stress  Waves  in  Laminated  Plates:  Using  similar  boundary 
conditions,  the  transient  stress  on  the  axis  of  symmetry  was  solved  during 
the  first  two  reflections  for  several  interesting  combinations  of  laminated 
plates. 

The  vertical  and  the  radial  stresses  as  a  function  of  time  and  position 

on  the  axis  of  symmetry  are  plotted  in  Figures  B-8  through  B-13  for  the 

following  cases.  The  two  plate  materials  are  steel  and  titanium.  The  total 

thickness  in  each  case  is  the  same  (—  =  0.8).  The  ratios  of  the  thicknesses 

a 

of  the  materials  vary  from  3-1  to  2:2  to  1:3.  The  titanium  plate  is  first 
exposed  as  the  upper  surface  (toward  the  bullet)  and  then  the  order  is 
reversed. 

The  use  of  a  laminated  plate  system  shows  (l)  an  impedance  match  be¬ 
tween  the  front  plate  and  the  projectile  increases  the  stress  generated  at 


$0 


a? 


Figure  B-9  STRESSES  T  •  A  LAMINATED  PLATE 
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Figure  B-li  STRESSE 
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the  impact  interface.  (2)  the  magnitudes  of  the  stresses  near  the  first  wave 
front  in  the  hack  plate  are  smaller  than  those  in  the  single  plate,  and 
(3)  higher  magnitudes  of  tensile  stresses  are  generated  in  the  second  and 
third  reflected  wave  fronts  passing  through  the  back  plate. 

B . 2 . 3  Inelastic  Stress  Wave  Propagation 

The  second  problem  solved  was  the  propagation  of  a  stress  wave  in  a  plate 
including  the  effects  of  plastic  flow.  The  stresses  generated  at  impact  are  a 
monotonically  increasing  function  of  the  impact  velocity,  for  impact  velocities 
larger  than  some  transition  velocity,  the  dominant  stress  in  the  plate  in  that 
impact  process  exceeds  the  static  yield  stress  and  is  in  the  plastic  range. 
Hopkins  and  Kolsky  (i960)  computed  the  transition  velocity  from  elastic  to 
plastic  behavior  for  aluminum- to- aluminum  impact  and  found  it  to  be  in  the 
neighborhood  of  40  feet,  per  second.  The  striking  velocities  of  interest  for 
the  ballistic  impact  studies  lie  above  this  value.  Therefore,  an  investigation 
into  stress  wave  propagation  in  the  plate  where  the  stresses  exceed  the  yield 
stress  of  the  material  was  conducted. 

A  study  of  plastic  flow  under  dynamic  and  transient  loading  conditions 
has  been  made  which  far  exceeded  the  yield  point  (Fugeiso  1962).  Using  the 
continuum  theory  of  microscopic  Imperfections  of  intra-crystalline  and  inter- 
crystalline  nature,  a  constitutive  stress-strain-strain  rate  equation  was 
developed.  This  constitutive  equation  well  represented  the  actual  material 
behavior  for  face-centered  and  body-centered  metals.  The  basic  strain  rate 
equation  given  by  this  theory  is 
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where  e . .  is  a  strain 

a. .  is  a  stress 
10 

is  a  function  of  the  stresses  and  temperature 
The  form  of  this  stress-strain-strain  rate  equation  is  equivalent  to  that  of  a 
visco-plastic  body.  The  actual  coefficients  of  the  body,  however,  are  nonlinear, 
as  indicated  by  the  theory  for  dislocation  movement.  Temperature  effects  are 
also  included  in  the  dislocation  theory. 

To  demonstrate  that  the  dislocation  model  and  its  predictions  indeed  give 
reasonable  stress-strain  curves,  the  stress-strain  curves  generated  at  constant 
strain  rate  under  various  boundary  conditions  are  shown.  The  stress  system  is 
biaxial.  In  both  problems  demonstrated,  the  longitudinal  stress  is  plotted  as 
a  function  of  the  longitudinal  strain.  The  longitudinal  strain  is  applied  at  a 
constant  rate  which  is  indicated  (Figure  B-l4).  In  the  first  case,  the  de¬ 
formation  corresponds  to  the  standard  tensile  test,  namely,  the  lateral  stress 
is  zero  and  the  body  is  free  to  deform  laterally.  In  this  case  the  stress 
starts  off  at  zero,  rises  linearly  until  a  yield  stress  is  reached,  at  which 
point  the  stress-strain  curve  bends  over  and  continues  almost  flat.  For 
other  models  of  microscopic  behavior,  a  reduction  in  the  stress  reached  at 
this  point  was.  obtained  (Figure  B-15). 

The  second  case  shown  is  the  longitudinal  stress  vs.  longitudinal  strain 
for  confined  deformation.  In  this  case,  the  lateral  strain  is  zero,  the  lateral 
stress  being  sufficient  to  confine  it  from  lateral  motion.  Again,  the  longi¬ 
tudinal  strain  is  applied  at  a  constant  rate.  Initially,  the  stress  varies 
linearly  with  the  strain,  the  modulus  being  equal  to  X  +  2j u.  When  the  yield 
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Figure  B-l4  STRESS  VS.  STRAIN  AS  A  FUNCTION  OF  STRAIN  RATE  3Y  TIE  THEORY  OF 
DISLOCATION  MOTION  -  CONSTANT  NUMBER  OF  DISLOCATIONS 
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Figure  B-15  STRESS  VS.  STRAIN  AS  A  FUNCTION  OF  STRAIN  RATE  BY  THE  THEORY 

OF  DISLOCATION  MOTION  -  FRANK-READ  RECENERATION  OF  DISLOCATIONS 


conditions  are  met,  the  longitudinal  stress  does  not  flatten  out,  but  now 
proceeds  as  a  function  of  the  longitudinal  strain,  increasing  at  a  rate  equal 
to  the  bulk  modulus  of  the  material.  The  lateral  stress  increases  at  the  same 
rate  after  yielding  and  has  the  same  sign  as  the  longitudinal  stress  (Figure  B-l6). 

In  these  cases, the  elastic  portion  of  the  deformation  is  linear.  The 
various  dislocation  models  were  chosen  to  represent  experimentally  observed  dis¬ 
location  behavior  in  solids.  In  Figure  B-1'7,  they  show  the  behavior  of  the  yield 
stress  (defined  as  the  limiting  stress)  for  the  unconfined  stress-strain  de¬ 
formation  at  constant  strain  rate  plotted  vs  the  strain  rate.  The  model  here  is 
that  of  titanium  alloy  yielding  by  dislocation  movement  given  by  the  Frank-Read 
source.  The  experimental  points  are  the  measured  values  of  Wentz  and  Catlin 
(1955)-  The  agreement  is  fairly  good. 

The  full  solution  of  the  equations  of  motion  for  thic  impact  problem  in¬ 
volves  the  computation  of  both  elastic  and  inelastic  behavior  in  at  least  two 
dimensions  and  time.  The  linear  elastic  solution  for  the  problem  is  already 
complicated  and  yet  it  is  the  only  model  which  will  allow  us  to  establish  an 
analytic  solution  (but  not  in  closed  form)  in  this  particular  geometry. 

Since  the  stresses  are  above  the  static  yield  point  for  most  of  the  deformation 
processes  involved  here,  an  approach  was  desired  which  demonstrated  the  effects 
of  the  inelastic  behavior  on  the  ballistic  resistance  and  on  the  properties  of 
the  perforation  and  penetration  process  m  general. 

Since  the  problem  is  the  impact  of  thin  plates,  thin  plate  approximation 
was  used  and  reduced  the  geometrical  problem  to  one  less  dimension  and  time; 
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Figure  B-1T  DYNAMIC  YIELD  STRESS  AS  A  FUNCTION  OF  STRAIN  RATE 


integrating  the  el  ‘ects  over  the  thickness.  Note,  from  the  elastic  solution, 
that  after  the  stress  wave  has  reflected  at  least  once,  the  vertical  variation 
of  the  stresses  over  the  plate  thickness  is  very  small.  Also,  the  magnitude  of 
the  vertical  stress  is  very  small. 

Proceeding  along  this  line,  a  mathematical  model  for  the  propagation  of 
transverse  deflections  of  a  thin  plate  including  plastic  flow  was  developed. 
Figure  B-l8  shows  the  projectile-plate  configuration.  The  only  variation  of  the 
displacement  allowed  was  in  the  vertical  direction,  an  approximation  which  is  a 
reasonable  description  of  the  impact  deformation  of  the  plate.  The  only  spatial 
variation  of  the  displacement,  strains,  and  particle  velocities  allowed  was  in 
the  radial  direction.  No  variations  in  the  vertical  direction  were  allowed. 

A  circular  plug  directly  underneath  the  cylindrical  projectile  was  assumed  to 
be  moving  at  a  velocity  equal  to  the  projectile  at  the  time  t  =  0.  This  accounted 
for  the  first  stress  wave  traverse  and  the  momentum  transfer  during  this  first 
instant  of  impact.  The  resulting  equations  and  boundary  conditions  are  : 

Momentum  equation 
dv  hr  t 

P  3t  =  57  +  7 


Constitutive  equation 


where 


P 

T 

7 

u 


0  (t ) 

0  (r) 


=  density  of  plate 

=  shear  stress 

=  shear  strain 

=  vertical  displacement 

=  vertical  velocity 

=  plastic  strain  rate  component 

=  A  e  form  chosen  to  match  experimental  data 


The  initial  conditions  and  boundary  conditions  are: 


where 


U  =  V  = 

0;  t  =  0; 

r  1 

R 

V  = 

vQ;  t  =  0; 

All 

R 

=  t  27rRh; 

r  ss 

R;  t  ^  0 

M  =  mass 

of  projectile 

m  =  mass  of  plug 
R  =  radius  of  projectile 
h  =  thickness  of  plate 


These  equations  of  motion  were  solved  numerically  on  MRD's  IBM  1620 
digital  computer.  A  finite  difference  method  for  the  solution  of  these 
partial  differential  equations  was  established,  checked  out  and  solved. 
Typical  results  for  the  velocity  and  the  shear  stress  as  a  function  of 
radial  position  and  time  are  shown  in  Figures  B-19  and  B-20.  In  this  par¬ 
ticular  case, an  initial  velocity  was  chosen  such  that  the  initial  shear 
stress  is  12  times  that  of  the  yield  stress  in  shear.  It  is  noted  that  the 
stress  at  the  wave  front  decays  more  rapidly  than  the  corresponding  ] inear 
elastic  solution  to  this  problem. 


?6 


5.00 


97 


The  solution  of  another  specific  problem  contains  information  of  interest 
to  this  particular  problem.  The  problem  solved  was  the  propagation  of  a 
spherical  wave  into  a  plastic  deforming  medium  (Fugelso  1962).  The  stress 
applied  at  the  boundary  of  a  spherical  cavity  was  above  the  yield  stress. 

It  was  found  that  the  stress  at  the  wave  front  which  elastically  decay  s~Tnver sely 
with  the  distance  from  the  cavity  decayed  at  a  far  greater  rate.  The  bulk  of 
the  disturbance  propagated  out  as  a  slower  velocity  and  a  large  permanent  de¬ 
formation  was  found  near  the  cavity  itself. 

B.3  Failure  of  the  Plate  and  Calculation  of  the  Ballistic  Coefficients 
B.3.I  Modes  of  Failure 

The  evaluation  of  the  stresses  in  the  plate  as  a  function  of  position  and 
time  is  an  important  step  in  the  evaluation  of  ballistic  performance.  If  the 
plate  is  to  fail,  it  must  fail  in  response  to  one  or  more  of  the  stresses. 

The  retardation  of  the  projectile  on  its  course  is  effected  by  the  stresses 
generated  at  the  contact  surface.  Therefore,  two  effects  in  determination  of 
resistance  to  perforation  must  be  established. 

1.  Does  the  plate  break?  If  so,  where  and  under  what  conditions? 

2.  If  the  plate  breaks,  how  much  is  the  projectile  retarded? 

It  is  conceivable  that  the  plate  could  break  and  the  bullet  would  not  per¬ 
forate  the  plate. 

The  description  of  the  retarding  and  stopping  of  the  projectile  by  the 
plate  was  made  by  integrating  the  stresses  that  act  over  the  surface  of  the 
projectile  into  a  resisting  force  acting  against  bhe  projectile.  The  model 
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of  stress  wave  propagation  and  failure  of  the  plate  material  and  of  the  pro¬ 
jectile  material,  together  with  the  geometry  of  the  plate  and  projectile,  will 
determine  the  magnitude  of  the  resisting  force  and  the  duration  of  the  resisting 
force. 

In  response  to  the  stresses  that  are  generated  in  the  plate  system  during 
the  impact,  the  plate  may  fail  "by  several  mechanisms.  The  mechanisms  are: 

1.  Fracture  in  the  Direct  Wave  -  the  plate  may  fracture  in 
response  to  high  stress  in  the  stress  spike  just  behind  the  first  stress 
wave  front. 

2.  Spall  -  the  plate  may  fracture  in  response  to  high  tensile 
vertical  stress  generated  by  the  reflection  of  the  first  stress  wave  from 
the  back  surface. 

3.  Radial  Fracture  -  the  plate  may  fracture  in  response  to  the 
tensile  stress  in  the  radial  direction  which  is  established  after  the 
decay  of  the  initial  stress  peak  at  the  wave  front. 

4.  Petalling  -  after  the  first  traverse  of  the  first  stress  wave 
the  plate  may  fracture  in  response  to  high  radial  and  circumferential 
tensile  stresses  that  are  developed. 

5.  Plugging  -  after  the  first  traverse  of  the  first  stress  wave 
the  plate  may  fracture  in  response  to  the  high  shear  stresses  developed 
in  the  plate  near  the  projectile. 

B.3.2  Calculation  of  the  Critical  Perforation  Velocity  and  the 
Residual  Velocity 

To  generate  a  prediction  for  ballistic  resistance  of  a  plate  to 
penetration,  examine  the  method  of  failure  in  the  direct  wave.  Consider  the 
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following  mechanism,  representative  of  the  direct  stress  wave  failure. 

Assume  that  the  plate  fails  in  compression,.  The  stress  at  the  wave  front 
directly  underneath  the  plate  in  the  vertical  direction  is  compressive. 

The  two  lateral  stresses,  namely,  the  radial  and  the  tangential  stress  on 
the  axis  of  symmetry,  are  equal  and  are  both  equal  to  a  given  fraction  of  the 
vertical  stress. 

Since  the  standard  material  test  for  failure  in  compression  has  no  lateral 
confining  pressure,  the  following  identification  between  the  standard  tensile 
failure  test  and  failure  is  made  here.  The  material  is  considered  to  have 
failed  if  the  maximum  shear  stresses  developed  are  equal  to  the  critical  shear 
stress.  In  this  case  the  difference  between  the  vertical  stress  and  the  radial 
stress,  which  is  equal  to  the  tangential  stress,  is  equal  then  to  the  tensile 
yield  strength  of  the  material.  (Figure  B- 21 ). 

At  the  back  side  of  the  plate,  the  vertical  stress  has  the  lowest  value 
of  the  compressive  stress  in  the  direct  wave.  The  compressive  stress  there 
is  a  function  of  the  following  three  parameters: 

1.  Impact  pressure  at  the  interface  between  the  bullet  and  the 
projectile.  The  impact  pressure  is  proportional  to  the  acoustic  impedance 
of  the  plate  projectile  combination  and  to  the  impact  velocity. 

2.  The  geometrical  attenuation  in  the  elastic  wave  propagation 
across  the  plate. 

3.  The  attenuation  by  plastic  flow. 

An  impact  velocity  may  be  computed  such  that  the  lowest  stress  in  the 
back  of  the  plate  is  exactly  equal  to  the  stress  required  for  fracture . 
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A  critical  velocity  for  perforation  is  thus  generated: 
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This  critical  velocity  is  a  function  of  the  yield  strength  of  the  plate,  the 
plate  geometry  ^given  through  the  ratio  of  h,  the  thickness  of  the  plate,  to 
the  bullet  radius,  and  to  the  plate  and  projectile  material  properties  such 
as  the  elastic  constants  and  density.  At  zero  plate  thickness, the  critical 
velocity  is  finite  and  ai  very  large  plate  thicknesses  It  approaches  a 
linear  increase  with  h.  The  residual  bullet  velocity  versus  striking  velocity 
was  calculated. 
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The  experimental  work  of  Zhurkov  and  Sanfirova  (1958)  has  indicated  that, 
at  least  in  constant  stress  (tension), there  is  a  relationship  between  the  time 
to  fracture  and  the  level  of  stress.  The  time  to  fracture  is  essentially  the 
inverse  of  the  plastic  component  of  the  strain  rate  corresponding  to  that  stress 
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as  determined  from  the  dislocation  theory  of  plastic  flow.  The  actual  magnitudes 

of  the  time  to  fracture  that  were  measured  by  Zhurkov  are  very  short  of 

-8 

the  order  of  10  seconds  for  very  high  stresses.  Figure  B-22  shows  an  ex¬ 
perimental  time  to  fracture  versus  stress  level  obtained  by  Zhurkov  in  his 
studies  on  fracture. 

For  a  first  approximation, the  time  to  fracture  at  a  stress  above  the  static 
yield  stress  of  the  material  is  so  short  that  it  can  be  ignored  and  the  time  to 
fracture  below  this  yield  stress  is  so  long  that  the  material  will  remain  un¬ 
broken.  As  an  approximation,  this  turns  out  to  be  a  very  good  description  for 
the  fracturing  and  failure  criteria  that  would  follow  the  first  wave  across  the 
plate . 

For  thin  plates,  the  time  to  fracture  is  comparable  to  the  transit  time  of 
an  elastic  wave.  For  a  steel  plate  .01  inch  thick,  the  time  that  a  longi¬ 
tudinal  stress  takes  to  traverse  across  and  back  is 


„  2  X  0.01  m 

t  ~  'j.^1  . — 7 -  =  10  seconds 

200,000  m/sec 

This  time  is  ten  tin.es  longer  than  the  time  to  fracture.  This  time  delay  should 
be  taken  into  account.  For  thinner  plates  correspondingly  larger  stresses  at 
impact  are  needed  to  effect  fracture.  (This  fact,  by  itself,  does  not  imply 

that  thinner  plates  are  better,  ballistically  speaking).  The  effect  of  time 
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delay  to  fracture,  using  ta  =  10  sec,  does  not  significantly  alter  the  pre¬ 
dictions  of  the  critical  perforation  velocity  for  plates  with  thickness 


h  =  0.01  inches. 


The  stresses  in  the  spike  decay  extremely  rapidly  by  dissipation  in  plastic 
flow.  Explicit  expressions  for  the  decay  of  a  stress  wave  by  plastic  flow  have 
been  computed  for  the  propagation  of  a  spherical  wave  (Fugelso  1962)  and  a 
cylindrically  expanding  wave  (Fugelso,  Arentz  and  Davidson  1962).  For  impact 
stresses  far  above  the  yield  stress  for  the  distance  the  stress  wave  must 
propagate  in  the  present  problem,  the  stresses  should  decay  by  a  factor  of  1.5 
to  2  by  the  time  the  stress  wave  reaches  the  back  cf  the  plate.  Thus,  the 
critical  perforation  velocity  estimate  should  be  raised  by  a  factor  of  1.5  to  2. 

After  the  first  stress  wave  propagates  across  the  plate,  it  reflects. 

Since  the  rear  surface  of  the  plate  is  free,  the  normal  or  vertical  stress  at 
the  reflected  wave  front  is  tensile  and  the  radial  stress  in  the  reflected 
wave  front  is  also  tensile.  These  stresses  will  add  to  the  existing  stresses 
left  after  the  decay  of  the  pressure  spike  in  the  direct  stress  wave.  The 
residual  normal  or  vertical  left  after  the  passage  of  the  spike  is  compressive 
and  the  residual  radial  stresses  in  the  plate  after  the  passage  of  the  spike 
are  tensile.  These  stresses  add  to  form  a  composite  stress  system. 

Since  both  of  these  stresses,  the  vertical  and  the  radial j are  tensile, 
the  failure  of  the  plate  must  be  considered  when  either  one  of  them  exceeds 
the  tensile  yield  stress.  Consider  the  radial  stress.  Near  the  back  surface 
of  the  plate,  the  radial  stress  due  the  direct  wave  is 


The  radial  stress  in  the  reflected  spike  is 
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stress  at  the  back  surface  is 
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The  critical  perforation  velocity  by  fracture  due  to  radial  stresses  is 
given  by  the  impact  velocity, which  makes  the  radial  stress  equal  to  the 
tensile  yield  stress  at  the  rear  surface  of  the  plate.  Thus, 


v  =  Z&  1  r  gty  r  2(i-y)r  1 

1/2 f-1  +  f"1/2  .  (l-v)  t2viL;'Z  i 


The  vertical  stress  in  the  reflected  stress  wave  is  also  tensile.  At  some 
distance  sufficiently  far  from  the  rear  surface  so  that  the  stress  spike  in  the 
direct  wave  has  decayed  to  zero,  the  vertical  stress  in  the  reflected  stress 
spike  is  (^denotes  this  distance) 


a 
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At  this  same  point,  the  residual  stress  in  the  direct  stress  wave  is 


The  vertical  stress  at  z  =  £*  is 


Making  the  approximation  (*•  )  «  (~) 
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the  tensile  vertical  stress  is  aTjproximately 
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The  critical  perforation  velocity  for  failure  in  response  to  large  normal  or 
vertical  reflected  stresses  is  obtained  by  determi;  \ng  the  impact  velocity  at 
which  the  normal  stress  is  equal  to  the  tensile  yield  stress. 


V  =  fk  _ -i _  _ % _ 

c  Z  2  (fV2.L) 

This  model  of  failure  is  failure  by  spall  and  is  normally  characterized  by  a 
plane  failure  near  the  back  side  of  the  plate.  A  small  group  of  pieces  be¬ 
tween  this  plane  and  the  back  surface  may  separate  and  fly  off  the  back  side. 

One  of  the  features  in  the  equations  for  the  propagation  of  the  stresses 
is  the  magnitude  of  the  stress  generated  at  the  impact  force.  The  exact  de¬ 
termination  of  the  magnitude  of  this  stress  is  a  most  difficult  analytic 
problem.  For  a  first  estimate, the  uniaxial  impact  of  two  linear  elastic 
bodies  was  used  to  estimate  the  stress  over  the  contact  area.  For  the 
generation  of  the  stresses  that  are  propagated  into  the  plate,  a  linear 
elastic  model  was  used  to  compute  the  stresses  generated  at  the  impact  area. 


In  this  model, the  impact  stress  (J^.  is 

„  .  P°C°  v 

I  ^ c,  +P0C0  s 

where  P^Pq  are  densities  of  the  two  materials 

Ci,C0  are  the  longitudinal  sonic  velocities  of  the 
two  materials 

V  is  the  striking  velocity  of  the  impact 
s 

An  additional  factor  in  the  impact  process  was  checked.  The  elastic-plastic 
stress-strain  curves  in  uniaxial  compression  were  used  to  compute  the  impact 
stresses.  The  elastic-plastic  curves  are  linear-elastic  in  volume  compression 
and  linear-elastic  fully-plastic  in  shear.  The  maximum  shear  stress  is  one- 
half  the  tensile  yield  stress  of  the  material.  The  stress-strain  curve  is 
shown  in  Figure  B-23  and  the  impact  stress  vs  impact  velocity  curve  for  uni¬ 
axial  deformation  is  shown  in  Figure  B-24.  The  impact  stresses  were  lower 

hy  about  10$  for  steel- aluminum  2024-T6  impact  at  V  =  4000  ft/sec.  The  in- 

s 

crease  in  the  critical  perforation  velocity  in  the  previous  modevjas  about 
10$  at  that  impact  velocity. 

If  the  deformation  at  the  impact  is  treated  as  unconfined  compression, 
the  maximum  stress  that  can  be  generated  is  0c^,  the  compressive  ultimate 
stress. 

Then  cl.  =  a 

I  cu 

In  addition  to  the  static  forces,  a  dynamic  pressure  term  proportional 
to  P-lV2  must  also  be  added  to  account  for  the  inertia  of  the  plate  material. 
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Figure  B-24  IMPACT  STRESS  VS.  IMPACT  VELOCITY 
FOR  UNAXLA1  IMFACT 


For  the  uniaxial  elastic  model  of  impact,  this  term  is  quite  small  compared  to 
the  stress  in  the  pcv  term  for  impact  velocities  up  to  $000  ft/sec  for  steel- 
aluminum  impact.  For  the  plastic  unconfined  model,  however,  the  inertial  term 
is  quite  large  and  cannot  be  neglected. 

If  the  plate  has  not  failed  after  the  first  traverse  of  the  first  stress 
wave,  the  projectile  may  still  be  moving  with  a  high  velocity,  particularly  in 
those  impacts  where  the  projectile  length  is  very  much  larger  than  the  plate 
thickness.  The  failure  of  the  plate  occurs  as  fracture  in  response  either  to 
the  large  tensile  membrane  stresses  developed  in  the  plate  underneath  the 
projectile  or  to  the  large  shear  stress  developed  in  the  plate  near  the  edge 
of  the  projectile. 

Nishiwaki  (1951 )  derived  the  equations  of  motion  for  the  normal  impact 
of  a  cylindrical  bullet  with  a  conical  cap  impacting  a  plate.  In  this  case, 
the  bullet  is  treated  as  a  rigid  body  moving  in  the  plate  with  some  initial 
velocity.  The  equations  of  motion  of  the  projectile  are  the  equations  of 
motion  of  a  rigid  body  acted  upon  by  time  and  position  dependent  force  and 
torques.  The  projectile  is  acted  upon  with  a  resisting  force  by  the  plate 
consisting  of  a  static  stress  given  by  the  particular  mode  of  failure  of  the 
plate  and  a  dynamic  stress  proportional  to  the  square  of  the  velocity.  The 
motion  of  the  projectile  is  determined  until  it  passes  through  the  plate  with¬ 
out  regard  to  the  details  of  the  deformation  of  the  plate.  In  this  time  of 
passage  the  forces  are  acting  upon  the  bullet  and  the  bullet  slows  down. 

The  perforation  parameters  of  ductile  metallic  plate  for  the  various 
models  of  plate  failure  are  calculated  by  using  different  forms  of  the  static 


component  of  the  resisting  force.  In  particular,  three  forms  of  the  static 
resisting  force  may  be  descriptive  of  the  plate  deformation.  For  very  thin 
plates,  the  stress  acting  against  the  projectile  may  be  a  function  of  the 
plate  thickness  as  for  the  shear  mode  computed  in  the  previous  section. 


J  =  a  (-) 
I  ty^a' 


The  resisting  stress  may  also  be  computed  as  a  function  of  plate  thickness  if 


the  membrane  stresses  in  the  plate  are  large: 


or  =  bo  (£) 

I  ty  a 

In  these  cases,  the  resisting  force  is  proportional  to  the  plate  thickness. 

For  moderately  thick  plates,  the  resisting  forces  of  the  plate  may  be  estimated 
as  the  compressive  ultimate  stress.  The  compressive  ultimate  stress  is  the 
largest  stress  that  can  be  developed  in  the  plate  since  the  material  is 
assumed  to  be  nearly  unconfined  and  flows  laterally  away  frcm  the  impact 
zone.  These  models  were  computed  for  the  prediction  of  the  perforation  of 


cylindrical  projectiles. 


Consider  an  axisymmetric  projectile  of  arbitrary  cross  section  im¬ 
pacting  a  plate  (Figure  B-25).  The  angle  of  impact  is  such  that  the  axis  of 
the  projectile  is  normal  to  the  plate.  The  projectile  has  density  p^,  mass 
M.  The  plate  has  density  p^,  thickness  h.  The  projectile  has  velocity  V 

normal  to  the  plate.  The  initial  velocity  is  V  .  In  this  analysis,  the 

s 

lateral  velocity  is  zero. 


Just  after  the  impact,  the  projectile  is  resisted  by  the  stresses  de¬ 
veloped  in  the  plate  at  the  contact  area  between  projectile  and  plate. 

The  material  in  the  plate  next  to  the  projectile  is  assumed  to  move  laterally 
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Figure  B-25  CONFIGURATION  OF  PLATE  AND  PROJECTILE 
FOR  NISHIWAKI  MODEL 
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out  of  the  projectile's  way  during  the  penetration.  Consider  a  differential 
element  of  the  projectile  surface,  dA  (Figure  B-26).  The  normal  to  this 
surface  makes  an  angle  a  with  the  horizontal.  The  resisting  stress  on  this 
element  dA  is  given  by  a  static  component  of  stress,  plus  the  hydrodynamic 
pressure. 

a  =  +P2}/2  Sin2  a 

The  resisting  forces  acu  in  a  direction  normal  to  the  surface  element. 
Since  the  projectile  is  axisymmetric,  the  force  component  resisting  motion 
of  the  projectile  can  be  integrated  over  dA.  The  horizontal  component  of  the 
force  vanishes  and  the  vertical  component  is 

dF  =  -  (aj  +  p2  Vr2  Sin2a)  Sin  a  dA 

where  Vn  is  the  instantaneous  velocity  of  the  projectile  at  the  surface. 

The  equation  of  motion  of  the  center  of  mass  of  the  projectile  is 

V&Z  =  -  f  (o  +  p  V  2  Sin2a)  Sina  dA 

dt  J  cu  2  n 

The  expression  for  the  deceleration  of  the  projectile  is  quite  general. 

It  will  be  solved  explicitly  under  the  following  assumptions: 

1.  The  projectile  does  not  deform  during  the  impact. 

This  implies  that  dA  will  be  known  from  the  shape  of  the 

projectile  and  the  amount  of  penetration.  It  then  also 

follows  that  V  =  V. 

n 


i  i  «: 


2.  The  resistance  of  the  plate  is  effective  only  as  long  as 

some  portion  of  the  projectile  surface  is  attempting  to  move 
into  a  volume  occupied  by  the  plate.  Frictional  forces  or 
shear  tractions  between  the  plate  and  projectile  are  ignored. 


The  residual  velocity  for  the  cylindrical  projectile  is  made  by  assuming 
the  projectile  has  a  cylindrical  body  capped  by  a  conical  nose  whose  apex  angle 
is  a  and  taking  the  limit  as  a->  90°  ♦ 


The  residual  velocity  V  for  the  cylindrical  projectile  is 

V  2  .  c2  (e-yh-l)  +  V  2  e-yh 
r  '  '  s 


where 


and 


where  L  is  the  length 


2P2 

PlL 

of  the  cylindrical  projectile 


The  critical  velocity  for  perforation  is  the  value  of  Vg  where  =  0. 

Thus 

V  =  c  (e^-l)1/2 
s 

To  extend  the  analysis  to  laminated  or  composite  plates,  the  resisting 
force  parameters  are  taken  as  a  function  of  position  corresponding  to  the 
material  initially  there. 


The  laminated  plate  system  consists  of  two  plates  of  thickness,  h^,  h^; 
densities,  P^,  P2;  and  ultimate  stresses  or^,  0^.  The  index  1  refers  to  the 
front  material  and  the  index  2  refers  to  the  back  material.  Let  H  ■  h^  +  h^. 
Let  po  be  the  density  of  the  projectile. 


For  the  blunt  nosed  cylinder,  which  is  the  limit  of  a  cone-shaped  cylinder 
with  apex  angle  equal  to  a  straight  angle,  the  equations  of  motion  of  the 
projectile  are 

MS  =  A  (ffi  +  Pi  y2)  1  “  L>2 

where  H  is  the  mass  of  the  projectile 

A  is  the  cross  section  area  of  the  projectile 
cLz 

V  =  ~  is  the  velocity  of  the  projectile 
and  z  is  depth  of  penetration  of  the  nose 


The  density  and  static  stress  on  the  right  hand  side  of  the  equation  of 
motion  are  given  by 


ai 

=  0 

Pi 

=  0 

z  <  0 

ai 

'A 

Pi 

=  pi 

0  <  z  <  h. 

ai 

II 

r\?  ; 

Pi 

p2 

hl  <  z  <  H 

a. 

i 

=  0 

Pi 

=  0 

z  >  H 

The  projectile  impacts  the  plate  with  velocity  V  .  The  upper  plate 

s 

retards  the  projectile  until  the  penetration  depth  is  h^.  Just  at  this  point, 
the  projectile  velocity  is  V  The  second  plate  then  retards  the  projectile 
until  the  penetration  distance,  z,  is  equal  to  H  =  h^  +  h^.  At  this  time,  the 
projectile  velocity  is 

The  projectile  velocity  is 


v2  -  C,2  (e-^.1)  +  vs2  e-Yl2 
V2  .  c*  (e-^-1)  ♦  Vri2  e'1*2 


0  <  z  <  hx 


h^  <  z  <  H 


with 


C  2 

al 

C1 

1 

*t>| 

h“ '  1 

c  ^ 

- 

C2 

P2 

v 

2Pl 

Y1 

=  PoL 

Y2 

2p2 

"  PoL 

Thus,  the  projectile  velocity  in  the  second  plate  is 


2  -Y  (z-hj 
C2  (e  d 


w  ?  -Y  (z-h  )  -Yh 

1  -1)  +  CL2e  *  1  (e  1  L-l) 


-Y  (z-h  )-Y  h 

+  V  e  d  1 
s 


hL  <  z  <  H 


The  critical  perforation  velocity  for  the  laminated  system  occurs  just 

when 


Vr  =  0  at  z  =  H  =  h^  +  hg 


Thus 


.  0  y.h.  Y.  (h.  -h_ ) 

V 2  =  C  2  (e  1  1  -  e  1  1  2  ) 


9  Y.h.  Y0h0 
+  C22(e  1  X(e  2  2-l)) 


The  analysis  of  the  ballistic  model  that  best  fits  the  data  may  be 
extended  to  the  case  of  oblique  impact.  This  has  been  done  for  oblique 
impact  of  bombs  into  soil  (Fugelso  and  Jepson  1963;  Fugelso,  Holmstrom, 
Bardy  and  Gut.teriberger  1965).  The  solution  to  the  oblique  problem  is 
necessarily  numerical. 


m 
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